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Dependence of Galaxy Quenching on Halo Mass and 
Distance from its Centre 
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ABSTRACT 

We study the dependence of star-formation quenching on galaxy mass and environ- 
ment, in the SDSS (z ~ 0.1) and the AEGIS (z ~ 1). It is crucial that we define 
quenching by low star-formation rate rather than by red colour, given that one third 
of the red galaxies are star forming. We address stellar mass M*, halo mass Mh, den- 
sity over the nearest N neighbours Sn, and distance to the halo centre D. The fraction 
of quenched galaxies appears more strongly correlated with Mh at fixed M* than with 
M* at fixed Mh, while for satellites quenching also depends on D. We present the 
M*-Mh relation for centrals at z ~ 1. At z ~ 1, the dependence of quenching on M* 
at fixed Mh is somewhat more pronouced than at z ~ 0, but the quenched fraction 
is low (10%) and the haloes are less massive. For satellites, M*-dependent quenching 
is noticeable at high D, suggesting a quenching dependence on sub-halo mass for re- 
cently captured satellites. At small D, where satellites likely fell in more than a few 
Gyr ago, quenching strongly depends on Mh, and not on M*. The Mh-dependence of 
quenching is consistent with theoretical wisdom where virial shock heating in massive 
haloes shuts down accretion and triggers ram-prcssurc stripping, causing quenching. 
The interpretation of S n is complicated by the fact that it depends on the number of 
observed group members compared to N, motivating the use of D as a better measure 
of local environment. 
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Galaxies can be crudely divided into two classes: blue 
star-forming disc galaxies and red ellipticals dom i nated 
by old stars (IStrateya et all 12001 iKauffmann et all 120031 : 
iBlanton et all 120031 ; iBaldrv et all 120041 ). The blue galax- 
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ies lie on a sequence in the star formation rate (SFR) 
stellar mass (M„) plane in that t he more massive galaxies 
form stars at a higher rate (e.g.. iBrinchmann et al.l 1200 
Salim et alj|2007t iNoeske et al.ll2007i ; iGilbank et aljboi " 
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since z ~ 1, but has decreased in normalisation (zero point) 
to the present day. This sequence is often called the "main 
sequence", but we will refer to it as the "sta r -form ing se- 
quence" (SF sequence) following ISalim et al.l (|2007l ) . In a 
plot of SFR vs. stellar mass, the intrinsically red galaxies 
tend to lie below the SF sequence. 

The galaxy bimodality is manifest not only in their 
star formation properties, but also in their masses and the 
number density of surrounding galaxies. In the local uni- 
verse, less dense environments usually feature less massive 
blue star-forming disc galaxies while denser surroundings 
like the cluster environment tend to host massive red and 
dead SO and elliptical galaxies (first meas ured by lOemlerl 
1 19741 ; iDavis fc Gelled Il97rl |Pressleilll980l) . Pure ellipticals 
tend to be the most massive of intrinsi cally red galaxies and 
resid e in the centres of large clusters (|Kormendv fc Bender! 
120121 ). 

Large surveys of hundreds of thousands of galaxies have 
made it possible to study density correlations with unprece- 
dented stati stical power. In t he Sloan Digital S ky Survey 
(SDSS - see I York et ai1l200Ch , iHogg et al.l (|2003l ) find that 
for the blue cloud in the colour-magnitude diagram, the 
environment density correlates with colour ( redder galax- 



study (jQuadri et al.l 120121 ') shows that this star- formation- 
density relation persists out to z ~ 1.8. 

Some have suggested that these rela t ions are driven 
by stellar mass (Ivan den Bosch et al.l |2008j; iMarcillac et al 



Noeske et al.l ( 20071 ) find that this sequence has persisted 200£ ; iGriitzbauch et al. 2011 Y while other s ( Cooper et al 



20ld at z ~ 1 and lPeng et aj]|2010l . 1201 ll at z ~ 0) show 



ies l i e in denser environments) (see als o Balogh et al.l 



20041; iBlanton et al.l |2 005al ; iBaldrv et all l2006h . Likewise 



Kauffmann et al.l (2004]) found that density correlates with 
several galaxy properties, and most strongly with mean 
SFR/M* (or specific star formation rate - SSFR), in that 
the average SFR/M, of galaxies is lower in denser en- 
vironments. Yet the SFR-density relation seems not to 
apply to the SF sequence itself but is manifested in 
the ratio of quenched galaxies to st ar forming galaxies 
dCarter et al.ll200ll ; iBalogh et all 12004 iTanaka et alj|2004 



Rines et al.l 12004 IWolf et all 120091 ; Ivon der Linden et~aH 
20ld : lPenget aljboid ). 

There is evidence that environment-related quenching 
operates differently between the most massive galaxies of 
a group or cluster, usually assumed to be residing at the 
centres of the group potential well (and hence are called 
centr als), and those th at are in orbit around them (satel- 
lites). |Hogg_et_alJ ((2003J find that for the red sequence in the 
colour-magnitude diagram of galaxies, the most and least lu- 
min ous galaxies, i.e., the centrals and satellites respectively 
- see lBerlind et al.l (|2005l ). lie in higher density environments 
while intermediately luminous g alaxies on the red sequence 
lie in a density " s addle ". Some (|van den Bosch et al.l 12008 ; 
IPeng et al.l 120101 . l201ll ) even suggest that density-related 
quenching affects only satellites and not centrals. 

Th e picture at z ~ 1 is similar. In the DEEP2 galaxy 
survey (|Davis et al.l 120031 ) , just as the c o lour bimodality 
persists out to z ~ 1 jLin et all Il999l; llm et all |200ll ; 
iBell et alj|2004 IWillmer et al.ll200alFaber et al.||2007l ). the 
mean colour-d ensity relation ([Cooper e t al. 2006) and the 
SSFR-density (|Cooper et al.ll2008l ) relation are also pr esent 
at z ~ 1, while, at these redshifts. iBundv et al.l |2006l ) find 
that red-sequence galaxies are still more abundant in dense 
environments than their star-forming counterparts. A recent 



that the SFR-density relation is separable from the colour- 
mass relation. It has also been a point of contention whether 
or not the "environment-related" quenching, which is usu- 
ally measured on group scales, is fundamentally driven by 
the pr esence of near neighbour s , or the halo mass or both 



(e.g.. | Weinmann et al 



2006a 



Skibba fc Shethl 120091 ; fskibbj 120091 : IWetzel et al.l |2012al : 



van den Bosch et al. I200S ; 
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20 111 ). Thus in order to understand the drivers 



of quenching, it is necessary to understand first how envi- 
ronment density, stellar mass and halo mass relate to each 
other, properly treating centrals and satellites separately in 
a coherent picture, as we attempt to do in this study. 



1.1 Density, stellar mass, halo mass and satellite 
distance 

To first order, haloes are self-similar in that they all have 
the same density profile p(r/R v i T ) in units of the virial radius 
regardless of their mass. Thus, if the density of a galaxy's en- 
vironment (i.e., the number density of other galaxies around 
it - let's call it S) follows the halo mass density p, the galaxy's 
environment density S is not expected to be correlated with 
its halo mass. However, on closer inspection, the 5-halo mass 
relation is more complex and depends strongly on the 5 mea- 
sure, the number of observed galaxies within a group, and 
the limiting mass/magnitude of the galaxies counted in a 
group. For example, environment density is often measured 
by the distance to a galaxy's JVth ne arest neighbour (for 
helpful reviews of de nsity measures see lCooper et al.ll200a ; 
iMuldrew et alj|2012j ). On reflection it is clear that such a 
density, Sn, has two modes depending on whether or not 
the galaxy is isolated or is embedded in a massive halo con- 
taining many satellites. For isolated galaxies (and for groups 
with fewer than N satellites), Sn is a measure of the prox- 
imity to the next big haloes, and is therefore a measure of 
field density ove r a large volume (this ha s been pointed out 
by for example, IWeinmann et al.ll2006al ). We call this the 
cross-halo mode of Sn- For a galaxy belonging to a massive 
halo, Sn measures the density of galaxies within its own 
halo. We call this the single-halo mode of Sn- In the single- 
halo mode the number density of galaxies is expected to 
correlate with the halo mass under the assumption that the 
distribution of satellites follows the dark matter distribu- 
tion w hich is described by something similar to the NFW 
profile |Navarro et al.lll997i ). For more massive haloes, the 
dark matter profile falls off less steeply with radius than for 
less massive haloes, so that if you are a central galaxy in a 
more massive halo, the probability of finding the iVth near- 
est satellite close to you is higher. Similarly, if you are a 
satellite in a large halo, the probability that you lie closer 
to the centre (and thus to other satellites) is higher in more 
massive haloes. 

Because of these dual modes of density, the 5jv-M* and 
Sn-M^ relations, and hence quenching in these spaces, will 
depend on the group membership of the galaxy. It turns out 
that most centrals reside in groups with fewer members than 
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the typical iV (3-5), so the interpretation of Sn for centrals 
is relatively straightforward. However, a significant number 
of satellites will reside in groups of more than N members 
and thus the interpretation of 5n will be not be simple. This 
motivates the search for a better measure of "environment" 
for satellites. We propose in this study the use of the relative 
projected distance of the satellite from its group centre (D). 

These complex relations have not been clearly under- 
stood to date although so me have made major steps toward 
this end (see for example IWilman et al.l |2010| who studied 
how halo mass rela tes to environme nt measures on differ- 
ent scales, and also lHaas et al.ll2012i who studied the rela- 
tion between several environment measures and halo mass 
in dark matter si mulations). Recent group catalogues by 
iGerke et ail (120051 ) (for the All- Waveleng th Extended Groth 
Strip International Survey - AEGIS) and lYang et al.1 <|2007h 
(for SDSS) make it possible to disentangle these relations 
and inform our intuition. 

Also important to quenching studies is the understand- 
ing of the relation be t ween a galaxy's stellar mass and 
halo mass. lYang et al.l 1)20081 ) used their group catalogue 
of the SDSS to show that stellar mass for centrals cor- 
relates strongly with halo mass up to a certain M» after 
which centrals of the similar massive M* can reside in a 
wide range of halo masses. This spread in halo masses at 
high M* presents an opportunity to discriminate between 
halo-dependent quenching and M*-dependent quenching. 

With a proper understanding of these relations between 
density, mass, and distance for centrals and satellites, we 
can distinguish between different mechanisms of quenching, 
which we will describe next. 



1.2 Quenching 

Several broad categories of quenching mechanisms have been 
proposed, with the mechanisms acting on central galaxies 
being different from the mechanisms acting on satellites. 

Quenching mechanisms for central galaxies fall into two 
classes, those that remove existing gas from central galax- 
ies, which we term "internal", and those that prevent new 
gas from falling in and cooling, which we term "exter- 
nal". Internal m echanisms include active galactic nucleu s 
(AGN) feedback (|Granato et alj|2004 ISpringel et al.ll2005l ), 
quasar-driven winds generated by collisions of galaxies 
(jSpringel et al.l 120051 ) , and stellar (radiation, winds and /or 
supernova) feedback (|Dekel fc Silkli 19861 ; iFall et al.ll2010h . 

Among the external mechanisms is what we call 
"halo quenching". In this picture, haloes more mas- 
sive than the critical mass for shock heating, M cr i t ~ 
few x 10 12 M(7) are capable of sustaining a virial 



shock (e.g., IBirnboim fc Dekell 120031 : iKeres et all 120051 ; 
iDekel fc Birnbo im 2006) . The heat ing by such an expanding 
virial shock (IBirnboim et al. 2007]) , by the energy of grav- 
itational infall (|Dekel fc Birnboimll2008l ). or by AGN feed- 
back that naturally couples to the hot dilute gas when it 
is present, shuts down the cold gas supply and suppresses 
star formation in galaxie s that reside in that massive halo 
(|Dekel fc Birnboimll2006h . For a population of haloes, this 
should not be interpreted as a sharp threshold but rather as 
a mass range extending over one or two decades where there 
is gradual growth of the fraction of quenched galaxies as a 
function of halo mass. 



Broadly speaking, the above proposed external mecha- 
nism for quenching centrals should depend on halo proper- 
ties rather than properties of the galaxies themselves. The 
internal mechanisms, on the other hand, depend on condi- 
tions or properties of the galaxy itself, such as black-hole 
mass, merger state, or star-formation rate. This distinction 
is useful for testing the importance of quenching mechanisms 
as we attempt to do in this study. However we recognise 
that quenching could be a combi nation of external and in - 
ternal mechanisms. For example, IDekel fc Birnboim! (|2006D 
suggested that halo quenching is possibly connected to the 
"radio mode" of AGN feedback from centra l black holes 
which inhibits central cooling in hot haloes (see lCroton et al.l 
I2006J ), 

As for satellite galaxies, the above mechanisms may 
have operated on them while they were still centrals, and 
may continue to operate on them to some degree after they 
become satellites by falling into a more massive halo. How- 
ever, there are additional quenching mechanisms that are 
unique to satellites only. One such mechanism is falling into 
a more massive and hotter halo than the satellite's original 
halo, impeding gas coolin g onto that satellite, usually re- 
ferred to as strangulation (lLarson et al. 19801; Balogh et al.l 



|2000| ; IDekel fc Birnboiml 120061 ; IChoi et alj |2009|). Another 
mechanism is ram pressure stripping ( Gunn fc Gottl Il972l ; 
lAbadi et al.lll999T ) which occurs only when the satellite or- 
bital velocity within the larger halo greatly exceeds its own 
internal dynamical velocities. Ram pressure stripping is ef- 
ficient where gas densities are higher, i.e., in the centr al re- 
gions of hos t haloes. Tidal stripping (|Read et al.ll2006l ) and 
harassment (|Moore et al.ll 19961 . 1 199ST ) are two more quench- 
ing mechanisms unique to satellites. Their efficiencies de- 
pend on the mass ratio of the intera ction, orbital eccentric - 
ity, and high environment densities (|Villalobos et al.ll2012l ). 

Among the first to explore the imp ortance of exter- 
nal vs. internal quenching was a study by IWeinmann et al.l 
(|2006al ) who found that halo mass was a better predictor 
of quenching for centrals and satellites than galaxy luminos- 
ity, pointing to external mechanisms of quenching. However, 
since the range of luminosity is large for a given M*, any 
sharp trend of quenching with M, may have been diluted, 
leaving the question of internal vs. external quenching still 
unc ertain. 

iKimm et al.l (|2009l ) also studied quenching in relation 
to M* and M^ but could not decide which was more impor- 
tant for centrals. For satellites, they found that they were 
equally important, but did not address the importance of 
environment density or central proximity. 

Two recent studies of quenching (|Peng et al.l |2010| , 

120111 ) propose an empirical model of quenching in order to 
explain the observed correlations between quenching, envi- 
ronment density and stellar mass. They propose that in- 
ternal quenching mechanisms for centrals should correlate 
with SFR. Without choosing a specific mechanism, they 
suggested that such a mechanism would be consistent with 
any kind of stellar feedback or AGN feedback. Through the 
SFR- Af* relation of star forming galaxies (i.e., the SFR se- 
quence), a correlation between quenching and M* is ob- 
served. In addition to this internal SFR-dependent quench- 
ing, the observed correlation between quenching and den- 
sity is assumed to be satellite-specific quenching of the sort 
discussed above. This idea of two kinds of mechanisms for 
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quenching, one mass-dependent and one satellite-specific, 
has be en implemented in s emi-analytical models (see for ex- 
ample, ICroton et aDl2006h . 

We build on this picture by putting everything in the 
context of the halo model, which likewise predicts differ- 
ent behaviour of quenching for satellites and centrals. How- 
ever, a difference is that the halo quenching model predicts 
stronger dependence on Mh than with M*, for both cen- 
trals and satellites, and this different quenching dependence 
proves to be observable. In the process, we also choose a 
better indicator of environmental density, D in preference 
to 5m, which provides additional insight for satellites. In ad- 
dition, we aim to study these quenching relations at z ~ 1 
where the number density of haloes above Mh.crit is lower, 
so halo quenching is not expected to be as important since 
fewer galaxies reside in large enough halos for virial shock 
quenching. Finally, we improve on previous studies by using 
dust-corrected SFR instead of colour to separate quenched 
galaxies from star forming galaxies. It turns out that the 
choice of colour versus SFR to classify a quenched galaxy 
yields different quenching behaviour. 



1.3 The goals of this paper 

In this paper, our goals are twofold. First, we aim to form 
an intuitive and consistent understanding of the complex 
relations between 5n, M* and Mh for central galaxies with 
the addition of D for satellites using the group information 
available for the SDSS and DEEP2 surveys. We demonstrate 
that keys to understanding these relations include separat- 
ing the relations for centrals and satellites and the fact that 
5n behaves differently in cross-halo and single-halo modes. 

With this understanding, our second goal is to test the 
physically motivated hypothesis th at quenching scales with 
Mh as predicted by halo quenching (IBirnboim fc Dekell 2003; 
iKeres et all 120051 ; [Dekel fc Birnboiml 120061 ). but differently 
for satellites and centrals when viewed in the diagrams of 
M,, Mh, 5n and D. Along the way, we compare quenching 
as measured by dust-corrected SFR with quenching as mea- 
sured by colour to see how much of an effect dust has on 
quenching behaviour. 

Our detailed mappings of quenched fractions as a func- 
tion of M*, Mh, 5n and D provide a rich repertoire of 
data for future testing of models. Matching all of these 
relationships quantitatively is a severe challenge. Several 
current semi-analytical models have problems reproducing 
som e aspects of the environmental dependencie s of quench- 
ing ( Wcin mann et al.l l2006bl ; iKimm et al.l 120091 ) and there- 
fore quenching trends need to be observed more closely. This 
is possible because star formation rates are becoming quan- 
titatively accurate under a wide range of conditions. 

The paper is organised as follows: In [J2] we describe the 
data sets, particularly the M*, environment density mea- 
sures, and the group catalogues used to estimate M h , and 
SFR measures for the SDSS and the Extended Groth Strip 
(EGS) field of DEEP2 (taken from the All- Wavelength EGS 
International Survey - AEGIS). In 331 we show that using 
dust-corrected SFR gives a better definition of star forming 
and quenched galaxies than rest-frame U — B colour. In ijlj 
we present and discuss the relations between Mh, Mm, 5n 
and D for the SDSS. Our main results are presented in Sj5] 
where we study the quenched fraction as a function of these 



quantities for central and satellite galaxies in the SDSS, and 
compare our results to previous work. In Sj6]we do the same 
exercise for the DEEP2 survey examining quenching in a 
context where halo quenching is not expected to be impor- 
tant. In SjZ] we summarise our results. 

All SFR and M* values are calibrated using the lKroupal 

(|200ll ) initial mass function (IMF), and concordance cos- 
mology (Ho = 70 km s _1 Mpc _1 ,f2M = 0.3, A = 0.7). All 
magnitudes are given in the AB system. Halo masses are 
converted to this cosmology with as = 0.9, fit = 0.04. 



2 THE DATA 
2.1 SDSS 

2.1.1 The sample 

The SDSS sample used throughout this analysis is limited to 
the DR6 area due to the coverage of the environment den- 
sity catalogue (to be described below). We limit the redshift 
range to 0.005 < z < 0.2 (abbreviated as < z < 0.2 or 
often loosely referred to as z ~ 0). After matching all cata- 
logues described below, the resulting sample size is 368,781. 
Matching all catalogues except the group catalogue, which is 
limited to M r — 5 log h = —19.5, the sample size is 459,174. 
This second sample is used in Sj3] but the first sample is used 
in the main analysis (S|4]and[5]|. 

We obtained ugriz photometry ("petro" values) 
and redshifts f rom t he NYU-VAGC (DR7) ca t alogue 
JBlanton et~aH l2005bl; I Adelman- McCarthy et al.l 120081 ; 
IPadmanabhan et al.l 120081 ) and used these data as i nput to 
the K-correction utilities of lBlanton fc Roweisl (|2007n (v4_2) 
to calculate Vniax from t he r-band limi t of the spectroscopic 
survey of 17.77 mag and iBesselll (|l990l ) U and B rest-frame 
absolute magnitudes (adjusted to h — 0.7). In the analyses 
to follow, we weight each galaxy by its 1/V max multiplied by 
the inverse of its spectroscopic completeness (also obtained 
from the NYU-VAGC). All quoted galaxy fractions and 
densities are weighted. 



2.1.2 SFR 

SFR esti mates are an upd at ed ve rsion of those de- 
rived in iBrinchmann et al.l (|2004f ). These new es- 
timates are taken from the DR7 catalogue and 
are provided online by J. Brinchmann et al. 
( http://www.mpa-garching.mpg.de/SDSS/DR7/ ) . As 

in Brinchmann ct al. (2004), these SFR estimates, cali- 
brated to the iKroupal (|200ll ) IMF, measure SFR within 
the fiber directly from Ha and H/3 lines for star forming 
galaxies. For galaxies with weak lines or showing evidence 
for AGN, they use the D4000 break calibrated to SFR using 
the star forming galaxies. Outsid e the fiber, the method fo r 
estimating SFR is improv ed over Brinchma nn et al.l i| 20041) . 
Following the method of ISalim et al.l ( 20071 ) . but without 
using UV flux points, they fit stochastic models of stellar 
populations to the observed photometry outside the fiber 
and constructed a probability distribution function (PDF) 
of SFR estimates. The means of the PDFs were added 
to the SFRs in the fiber for a total estimate of SFR. J. 
Brinchmann (private communication) estimate the typical 
error of these SFRs to be about 0.4 dex for star forming 
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galaxies, 0.7 dex for more quiescent galaxies and growing 
to 1 dex or more for dead galaxies (but for dead galaxies, 
these SFRs are more upper limits than measurements). Our 
own analysi s comp aring these SFR estimates with those of 
ISalim et all |2007f ) . which are derived from SED fitting of 
UV and optical light from GALEX, indicate that the error 
is closer to 0.2 dex for star forming galaxies. 

Our results are not changed s ignifi cantly if we use the 
SFR estimates from ISalim et al.l (|2007l ) since our analysis 
depends only on a division between star forming and passive 
galaxies ( jj3{ rather than on an absolute SFR value. However 
since the lSalim et al.l (|2007l ) SFR estimates cover only part 
of the DR4 (or DR7) a rea, we opted to use the updated 
iBrinchmann et al.l ((20041) SFR estimates which increase our 
sample size 10-fold. 



2.1.4 Stellar masses 

Stellar mass M* is taken from the M* cata- 
logue provided online by J. Brinchmann et al. 
(http://www.mpa-garching.mi3g.de/SDSS/DR7/ ). They 
derive M,, calibrated to the JKroupal (2001) IMF, t hrough 
SED fitting similar to the method that ISalim et al.l (|2007l ) 
used for estimating SFR's, but without using UV con- 
straints. The m et hod i s also related to that used by 
iKauffmann et al.l (J2003), who fitted spectral features 
inste ad of photometry. These stellar mass estimates differ 
from IKauffmann et al.l (|20031 ) by less than 0.1 dex for 
M» > 10 9 Mq for more than two-thirds of the galaxies. The 
formal 1-cr errors derived from the 95% confidence intervals 
of the probability distribution are typically about 0.05 dex 
or less. 



2.1.3 The group catalogue and halo masses 

lYang et al.l (|2007l ) constructed a group catalogue and calcu- 
lated group halo masses for the SDSS DR4 sample. Yang et. 
al. also appl ied the same algor ithm to the DR7 sample, as 
described in lYang et al.l (|2012l ) , and we use this catalogue 
throughout this study. We briefly describe their group finder 
and calculation of halo masses below. 

The group finder consists of an iterative procedure that 
starts with an initial guess of group centres and member- 
ship (the result of a friends-of-friends algorithm) . For each of 
these tentative groups they calculated the group size, mass 
and velocity dispersion and used these to calculate a pro- 
file for the number density contrast of dark matter part i- 
cles based on a spherical NF W profile (|Navarro et al.ll 19971 ) . 
Group membership was then recalculated based on the num- 
ber density contrast expected at the distance of the potential 
group member to the group centre. Then the centres, sizes, 
masses and velocity dispersions were recalculated using the 
new membership, applying completeness corrections to ac- 
count for missing members. Then the process was repeated 
until there were no further changes to group memberships. 
They tested their group finder on an SDSS mock catalogue 
found that their group finder successfully selected more than 

90% of true haloes more massive than 1O 12 M0. 

Once the group catalogue was constructed, lYang et al] 

(|2012f ) assigned halo masses using an abundance matching 
method. They rank-ordered the groups by group stellar mass 
and ass igned halo ma s ses ac cording to the halo mass func- 
tion of Warren et al.l (120060 with the transfer function of 
lEisenstein fc Hul ( 19980 . Comparing these halo masses with 
groups found in a mock catalogue, they find a rms scat- 
ter of about 0.3 dex. They also calculate halo masses by 
rank-ordering group luminosities, but we choose to use the 
rank-ordered group masses because stellar mass is a bet- 
ter indicator of h alo mass than luminosity (see for example 
iMore et al.1l201lh . 

Using this group catalogue we defined the most massive 
member in each group to be the central galaxy of the group 
and all other members to be satellites. Galaxies that were 
deemed by the group finder to be isolated are also defined 
as centrals. 



2.1.5 Environment density 

The environment de nsity measure log(l + 8n) is defined in 
ICooper et al.l (|2005l ). but we briefly describe it here. The 
number density (per unit area) of a galaxy is calculated from 
the projected distance to the iVth nearest neighbour (within 
a velocity window of ±1500 km s -1 to exclude foreground 
and background objects) and then divided by the median 
surface density of the whole sample at the galaxy's redshift 
to give an overdensity relative to the median density 1 + 5n ■ 
This median density is calculated in bins of Az = 0.02. Edge 
effects are minimised by removing all galaxies within 1ft" 1 
comoving Mp c from the survey's edges. For the SDSS (DR6) 
ICooper et al.l (|2009l ) calculated 1 + 8n for N — 5 and these 
are the density estimates used here. Typical 1-<t uncertain- 
ties for this density measure are about 0.5 dex. This environ- 
ment density measure was computed for the galaxies of the 
SDSS DR6 photometric catalog which is limited to objects 
brigh er than r — 22.2 in asinh magnitudes l|Lupton et al.l 
Il999h . 

In this work we also define the relative distance D 
of satellites to the centres of their haloes as a mea- 
sure of environment density. Using the group catalogue in 
32,1.31 to define centrals and satellites, we calculate the 
projected distance d pr oj of each satellite to the central 
galaxy of its group and divide that b y the virial radius 
flvir = 120(M h /10 11 M Q ) 1/3 kpc (e.g.. iDekel fc Birnboiml 
2006). Thus D = d ploi /R vir . 



2.2 AEGIS 

2.2.1 The sample 

AEGIS is a collaboration between many survey teams us- 
ing a variety of space- and ground-based facilities to ob- 
serve the Extended Groth Strip (EGS) field over var- 
ied areas (| Davis et al.l 120071 . Newman et al., in prepara- 
tion). Our sample is limited to areas cov ered simultaneousl y 
by Spitzer MIPS 24pm p hotometry (iDayis et al.l 120071 ') , 
CFHT BRI photom etry (ICoil et al.l 120041) . DEEP2 and 
DEEP3 spectroscopy (|Davis et al.ll2007l; I Cooper et alll201ll . 
l2012f ). Palomar A-band phot ometry dBundv et al.l 120061 ). 
and GALEX UV photometry (|Salim et al.l 120091 ). Further- 
more, only galaxies with redshift quality 3 or 4 (which means 
that at least two certain spectral features were identified by 
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eye - see lDavis et al.ll2007l ) were included this analysis. Ob- 
jects also had to have a probability of being a galaxy of 
greater than 0.4 (unresolved galaxies are assigned a prob- 
ability of 0-1 of being a galaxy base d on size, colour and 
magnitude - see iNewman et al.l I2012T ) . The sample size for 
these overlapping data sets for our redshift range of interest 
(0.75 < z < 1) is 1846. 

Using the Tf-correction utilities of iBlanton fc Roweia 

J2007J) (v4_2) we calculated V max from the 7?-band limit 
of the spectroscopic survey of 24.1 mag. We weight each 
AEGIS galaxy in this analysis by 1/Vmax multip li ed by 
the "optimal" weights calculated by IWillmer et al.l (| 20061 ) 
which account for co lour-dependent redshift success rate (see 
IWillmer et~alll2006l for d etails). 

Willmer et al.l (|2006j) also computed U and B absolute 



rest-frame magnitudes using their own Tf-corrections, and 
these are the magnitudes used in this analysis, adjusted for 
h = 0.7. 



2.2.2 SFR 

Several star-formation tr acers are available from the multi- 
wavelength AEGIS field (|Davis et al.ll2007h : 

• IR light from reprocessed UV photons from young, mas- 
sive stars that are dust-enshrouded. The rest-frame 15/im lu- 
minosi ty of galaxies is tightly correlated with total IR lumi- 
nosity (he Floc'h et alJl2005l : lBell et al.ll2005l . and references 
there in) which can be used to estimate SFR (jKennicutg 
1 19981 ). The IR data are from the Spitzer MIPS 24/xm pho- 
tometry (which corresponds to rest-frame 12-15/xm), which 
are used to derive the total IR luminosity using SEP tem- 
plates of I Chary fc Elbazl (|200ll ). lLe Floc'h et al.l J2005I ) find 



that the dispersion between IR luminosities estimated from 
several different SED libraries can reach up to 0.2 dex. 

• Non-extinction-corrected rest-frame UV luminosity 
from you ng, massive stars, derived from the observed B flux. 
Following [Bell et all (|2005l ). the rest-frame 2800A luminos- 
ity is used to estimate the total integrated 1216-3000A UV 
luminosity, which comes from ~ 100 Myr old populations. 
The UV luminosi ty is in turn used to estimate SFR following 
IKennicuttl (|l998l ). 

• Extinction-corrected UV luminosity derived from 
GALEX which observes in the near UV (NUV). Combined 
with the CFHT u* band, one may calculate the UV slope 
/3 between the rest-frame far UV (FUV) and NUV. Each 
galaxy is corrected for extinctio n using the correla tion be- 
tween /3 and FUV attenuation (|Seibert et al.ll2005l ). (How- 
ever, note that th e UV slope may also depend on age a nd not 
only on dust - see lKong et al.ll2004l ; ls"alim et al.ll2009l . We do 
not attempt to correct for this.) The FUV l umino sity is then 
used to calculate SFR following IKennicuttl (1 19981). This SFR 
deriv ation is described in more detail in ( Schiminovich et al.l 
120071 ). 

• Ha, H/3, [O II] emissio n lines from ionising photons from 
very massive, young stars (|Weiner et al.ll2007l ). For DEEP2 
spectroscopy, the Ha line is visible out to z ~ 0.37 (outside 
our redshift range of interest, 0.75 < z < 1), the H/3 line is 
visible at 0.35 < z < 0.86, and the [O II] line is visible at 
0.77 < z < 1.42. Only those galaxies with S/N > 2 are con- 
sidered a detected emission line SFR. (Note that line SFR's 
carry some limitations including the uncertainties of esti- 



mating dust extinction for the Balmer lines a nd the depen- 
dence of the oxygen lines on metallicity; see iMostek et al.l 
|2012i for details.) We used the mass-dependent flux ratios 
Fnp/Fna and Fon/Fna (for galaxies for which these lines 
observable) to estimate Ha from H,9 and [O II] at higher 
redshift s, and also to estimate the reddening correction 
Weiner et al.l 120071 ). Then the prescription of IKennicuttl 



1998) is applied to calculate SFR from Ha. 



ISalim et al.l (j2009f ) and Noeske et al. (in preparation) 
find that the combined SFR estimates from IR and non- 
extinction-corrected UV (IR+UV) yield good agreement 
on average with extinction-corrected SFR estimates from 
GALEX, and, where available, estimates from Ha. Emission 
line indicators yield SFR estimates that also agree with the 
other indicators except that they are systematically higher 
(0.2-0.4 dex). 

With this wealth of SFR indicators, we created a fidu- 
cial sample combining several indicators in the following way 
that we believe t akes advantag e of the strengths of each 
indicator (see also lNoeskell2012l ). Since IR+non-extinction- 
corrected UV (hereafter IR+UV for brevity) is growing in 
popularity as a reliable indicator of star formation (e.g. , 
lElbaz et al.1l2007l ; iNoeske et alll2007l ; ISalim et al]|2009l ) , we 
use IR+UV as our first estimate of SFR, and where this is 
not available, we use the extinction-corrected estimates from 
GALEX. Where neither IR+UV nor GALEX UV are avail- 
able, we use the SFR estimates from emission lines. In the 
redshift region where the ranges for detecting H/3 and [O II] 
overlap, we select H/3 as the preferred SFR indicator. Since 
these indicators produce SFR estimates that are offset from 
each other, we normalise the other indicators to IR+UV so 
that their average offsets are zero. To calculate these off- 
sets, we used a wide range in log SFR of -1.2 to 3. Overall, 
27% of this mix is derived from IR+UV, 43% from GALEX 
UV, and 30% from emission lines. We checked whether re- 
stricting the SFR sample to any single method used in the 
fiducial combination makes any difference. Results were not 
affected, and we do not discuss these other approaches any 
further. 



2.2.3 The group catalogue and halo masses 

iGerke et al.l (|2012j) constructed a group catalogue for the 
AEGIS field using a Voronoi-Delaun ay method (VDM) 
group finder (first implemented by iMarinoni fc Piazzal 
120031 ). This method makes use of a Voronoi partition of the 
galaxy catalogue into unique polyhedrons each containing 
one galaxy, and it also makes use of the Delaunay mesh 
which is a system of line segments connecting neighbour- 
ing galaxies. Starting with the smallest Voronoi volumes, 
the group finder finds all the neighbouring galaxies to a test 
galaxy via the Delaunay mesh within a certain cylinder. The 
size of the cylinder is determined iteratively as a function of 
the number of group members, and its free parameters are 
optimised to produce the bes t success in a mock AEGIS cat- 
alogue (described in detail in lGerke et al.l (|2012T l). Success is 
measured by completeness and purity defined as follows. One 
can check the mock catalogue to see if a group matches the 
halo that contains a plurality of the group's members and 
if a plurality of the halo's members are part of the group. 
Such a situation is considered a two-way match. Complete- 



Quenching by Halo Environment 7 



ness is the fraction of haloes that have a two-way match 
to a group found by the VDM group finder, and purity is 
the fraction of VDM groups that have a two-way match to 
a halo in the mock catalogue. The AEGIS group catalogue 
has a completeness of ~ 74% and a purity of ~ 67%. 

Using this group catalogue, we defined the most mas- 
sive member in each group to be the central galaxy of the 
group and all other members to be satellites. We also de- 
fined as centrals those galaxies that were deemed by the 
group finder to be isolated. To test the effectiveness of mas- 
siveness as an identifier of central and satellite galaxies, we 
identified the most massive members of the groups found 
by the group finder on 40 AEGIS mock catalogues. When 
compared to the mocks' underlying dark matter distribution 
(from the N-body simulation), we find that massiveness cor- 
rectly identifies a central or satellite on average 87±1% of 
the time, the ±1% being the standard deviation around this 
mean over the 40 mocks. The fraction of centrals identified 
in this way that are actually satellites is 12 ± 1% and the 
fraction of satellites that are actually centrals is 29 ± 1%. 

In order to estimate host halo masses, we assumed that 
each group roughly represents a distinct halo and estimated 
the halo masses of each group using an "abundance m atch- 
ing" method similar to that used bv lYang et al.l (|2007l ) (de- 
scribed in £|2 . 1 p - For each group, we calculated the total stel- 
lar mass M*,tot for all group members above M* = 10 M© 
in each group. Each group me mber was weighted by the 
"optimal" weights computed bv lWillmer et al.l (|2006l ) to ac- 
count for the spectroscopic completeness of the AEGIS sur- 
vey. Then we self-calibr ated the groups to account for miss- 
ing members following lYang et al.l (|2007l ). Specifically, we 
selected groups found in AEGIS between 0.2 < z < 0.4 and 
calculated the fraction of the group mass that is contributed 
by members above the stellar mass limits of (10,10.1,10.2) 
corresponding to t he limits at the redshifts (0.7,0.9,1.1) 
(jBundv et alj|200q ) . We then fit an exponential function to 
these fractions as a function of group mass for each redshift 
bin, and these exponentials became the redshift- and mass- 
dep endent fractions b y which we divided each group mass 
(see lYang et al.1120071 for further details). 

We then assigned a halo mass to each group by rank 
ordering M»,tot and matching them with the predicted dis- 
tribution of halo masses in redshift bins of 0.1. The halo 

using the mass 



mass function is that of ITinker et al 



transfer function of lEisenstein fc Hu 



), computed us- 
ing a code kindly provided by M. Cacciato. This method as- 
sumes that the halo mass function and the luminosity func- 
tion of groups match one-to-one. This is a rather simplistic 
assumption, but we are interested mainly in average trends 
with halo mass, and we expect that the relative masses be- 
tween groups of very different numbers of members to be 
quite robust. 

Using the same algorithm on the groups found by the 
group finder on the AEGIS mock catalogues, we find that 
the halo masses agree well with those of the mocks with a 
mean rms scatter of 0.37 dex for the rank-ordered stellar 



2.2.4 Stellar masses 

The stellar masses M* are derived in iBundv et al.l l|2006h . 
These authors use BRIK colours, as well as the J-band from 
Palomar when available in the EGS, and spectroscopic red- 
shifts to fit observed SED s to synthetic SEDs from the stel- 
lar population models of iBruzual fc Charlotl (J2003J). These 
model SEDs span a large range of star formation histories, 
ages, metallicities and dust content. From the fits, they con- 
struct a PDF for stellar mass-to-light ratio M*/Lk for each 
galaxy and use the median of the PDF as the best esti- 
mate mass-to-light ratio to calculate M*. Uncertainties in 
the model fitting combined with the uncertainties of the ob- 
served K-band luminosity yield typical final 1-cr uncertain- 
ties of 0.2-0.3 dex. 

In order to increase the statistical power of the AEGIS 
sample, we supplemented these stellar ma sses with those cal- 
culate d using the K-correction utilities of lBlanton fc Roweisl 
(J2007J) (v4_2) whenever Tf-band stellar masses were not 
available. Since the K-correction stellar masses were on av- 
erage smaller than the K-band masses by 0.05 dex we added 
this value to all the K-correction masses used in the sample. 
Supplementing the stellar masses in this manner adds 425 
galaxies to the 0.75 < z < 1 sample (for a total of 1846). 



2.2.5 Environment density 

ICooper et ail (|2008h calculated log(l + S N ) for the EGS 
field using TV = 3 and a velocity window of ±1000kms _1 , 
and these are the density estimates used here. Log(l + Sn) 
is described abo v e in |j2.1l Using the mock catalogues of 
lYan et al.1 (J2004l ). ICooper et all (|2005r i compared the effec- 
tiveness and accuracy of several measures of environment 
density and found that the projected distance to the third 
nearest neighbour provided the most accurate measure of 
environment density for the DEEP2 survey. Edge effects are 
minimised by removing all galaxies within l/i comoving 
Mpc from the survey's edges. Although the TV we use is dif- 
ferent for the SDSS and for AEGIS, the number densities 
used to calculate log(l + Sn) are higher in the SDSS than in 
AEGIS, and so the distances probed are comparable. The 
typical distance to the third nearest neighbour in AEGIS is 
about 1.5 h^ 1 comoving Mpc while the typical distance to 
the fifth nearest neighbour in SDSS is about 1.7 h^ 1 comov- 
ing Mpc. Typical 1-cr uncertainties for this density measure 
are about 0.5 dex. 



3 BIMODALITY IN SFR VERSUS COLOR 

The galaxy bimodality is usually discussed in terms of the bi- 
modality in the colour-magnitude diagram. Fig. [T] shows the 
colour-magnitude diagram of the SDSS sam ple with t he lin e 
dividing the colour bimodality used in lCooper et al.l |2008): 



U-B= -0.032 * (Mb + 21.62) + 1.175. 



(1) 



After matching all catalogs listed in £12.2.11 our sample 
includes 274 groups with more than one member, and 1361 
isolated galaxies. 



The red sequence, which lies above this line, is assumed to 
consist mostly of galaxies that have ceased to form stars 
while the blue cloud (below the line) consists of currently 
star-forming galaxies. Such a division is convenient because 
photometry of large samples is easy to do with relatively 
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Figure 1. The colour-magnitude diagram of the SDSS galaxies 
(0 < z < 0.2). The solid line divides the colour bim odality and is 
defin ed as U - B = -0.032 * (M B + 21.62) + 1.175 jCooper et alj 
120081) . The size of a pixel used to smooth the diagram is indicated 
at bottom left. 



good accuracy. However, techniques have improved for de- 
termining star formation rates for large samples, and we can 
now divide the galaxies by star formation directly rather 
than by colour as a proxy. 

We show the "physical" version (i.e., reflecting the phys- 
ical processes of star formation) of the colour-magnitude di- 
agram, namely the SFR-Af* diagram for the SDSS and the 
AEGIS in Fig. [2| The star-forming g alaxies lie on a sequence 
as described bv lNoeske et al.l (|2007T l . which we call the "SF 
sequence," while there is a population of galaxies lying un- 
der it, which we interchangeably call "low-SFR", "passive" 
or "quenched" galaxies. We defined the division between 
these populations in the following way. We assigned an ini- 
tial line to divide the population by eye. Then we performed 
a weighted linear least-squares fit to the points above that 
line and calculated the distance of each point from the line. 
After a Gaussian fit to these distances, we then exclude from 
the SF sequence those points that lie more than 2.5a below 
the line, and re-perform the fit. The process is repeated until 
no more points are excluded. Quenched galaxies are galaxies 
that do not lie in the SF sequence as defined in this way. We 
show the resulting divisions in Fig. [5] with the equation for 
the dividing line (solid black) being: 



log SFR = a log M, 



(2) 



where a = [0.64,0.31] and b = [-7.22,-3.04] for the redshift 
bins [0 < z < 0.2, 0.75 < z < 1]. 

When we refer to ratios of the low-SFR galaxies to the 
total number of galaxies, we call this the "quenched fraction" 
rather than the "red fraction" to remind the reader that we 
divide the galaxies by SFR and not by colour. 

In order to compare the SFR bimodality with the colour 
bimodality, we show in Fig.[3]the residuals from the SFR di- 
vision A log SFR = log SFR - Eq. ([2| versus the residuals 
from the colour division A(U — B) = (U — B)— Eq. JTJ for 
the SDSS. The horizontal solid line in Fig. [3] corresponds 
to the solid black line in Fig. [2}i. As explained above, we 



call everything above this line a star forming galaxy, or SF 
sequence galaxy, while everything below it is quenched. Sim- 
ilarly, the vertical solid line in Fig. [3] corresponds to the line 
in Fig. [l] that divides the red sequence (right of the vertical 
solid line in Fig. [3} from the blue cloud (left of the verti- 
cal solid line). The two horizontal dashed lines are ±0.15 
around the SFR division and are meant to guide the eye to 
the area of intermediate SFR. The vertical dashed lines are 
±0.05 around the colour division and mark the area of the 
green valley. 

From Fig. [3] we observe that the SFR bimodality mostly 
corresponds to the colour bimodality, but there are some no- 
table differences. In particular, green valley galaxies are not 
necessarily intermediate in SFR. Rather about 80% of green 
valley galaxies are still star forming (lying above the hori- 
zonal solid line). Furthermore, about 30% of red sequence 
galaxies lie in the SF sequence (these are star- forming and 
are pro bably reddened b y dust - see iBrammer et alj 12009 



and also lLotz et al.ll2008l who find a similar result for z ~ 1) 
whereas only 4% of the blue cloud lie below the SF se- 
quence (these may be post-starburst galaxies, those that 
exhibit no signs of ongoing star formation but experienced 
rec ent star formation with i n the last ~Gyr - se e for exam- 
ple , iDressler fc Gunnll 19831 ; iQuintero et~aJll2004 l . Therefore 
dividing the galaxies by U — B colour does not accurately di- 
vide the population by star formation mainly because dusty 
star formers appear red. We recognise that dust may not ex- 
plain fully explain the large population of star forming red 
galaxies since U — B and SFR contain large scatter and are 
constructed differently. However, the s e find ings are consis- 
tent with the results of iMaller et alj (J2009J), who find that 
about a third of red sequence galaxies in a volume-limited 
sample of the SDSS at 0.02 < z < 0.22 actually lie in the 
blue cloud after applying an inclination correction to their 
M g — M r colours, implying a dusty origin to the observed 
red colours of these galaxies. 

These dusty star formers are also the more massive 
members of the SFR sequence - the dusty star formers 
are on average 0.4 dex more massive than the entire SF 
sequence population. This is understandable as the more 
mass ive members of the SF sequence are more metal- 
rich l|Tremonti et alj 12004 ) and therefore have proportion- 
ately more dust. This is consistent with the results of 
IMaller et alj (|2009l ) who find that the inclination-dependent 
colour correction also depends on absolute K magnitude in 
that brighter galaxies have higher attenuation than dimmer 
galaxies of the same inclination angle. The fact that mas- 
sive star-forming galaxies are dusty will cause more massive 
galaxies to appear quenched if one defines quenching only 
by colour. This will be important in comparing our results 
with those of others who define quenching by colour only. 
We discuss this further in £15.31 



4 RELATIONS BETWEEN MASSES AND 
ENVIRONMENT 

In this section, we describe the complex relations between 
stellar mass M* , halo mass Mh , environment density Sn and 
satellite group-centric distance D for central and satellite 
galaxies in order to prepare the reader for a discussion of 
quenching as a function of these quantities in the follow- 
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Figure 2. Log SFR vs. log M* for galaxies in two redshift bins. The dashed blue line is a fit to the SF sequence. The solid black line 
shows the division between the "SF sequence" of galaxies (above the line) and the "low-SFR" or passive galaxies. The text describes 
how these lines are determined. The weighted ratio of low-SFR galaxies to the total number of galaxies is the "quenched fraction" . 



ing section. Summarising this section in brief, Mh correlates 
with M* for centrals at low masses, but shows a large spread 
at high M*, offering a chance to separate quenching effects 
between M 4 and Afh for centrals. Sn for centrals is a simple 
two-mode function of Mh while for satellites Sn is given by a 
complicated interplay between Mh and the location within 
a halo. The dividing line of interpretation of Sn is where 
the number of observed group members is less than N (the 
cross-halo mode) or greater than N (the single-halo mode). 
Due to the complex nature of 5m, the better measure of lo- 
cal "environment" for satellites may be radial distance from 
the group centre which is highly anti-correlated with Sn, 
but only in the single-halo mode. Readers interested mainly 
in the results of quenching as a function of M t , Mh and 
<5iv/radial distance may skip to Sj5] 

Note that in all the figures that follow for the SDSS 
analysis, one quantity is plotted as contours as a function 
of the x- and y-axes. These contours are computed in the 
following way: values of this third quantity are computed in 
pixels of the indicated size in the plots as long as the pixel 
contains more than m galaxies ("good" pixels), m is equal 
to 2% of the average number of galaxies in each pixel which 
contains at least one galaxy. The values are then smoothed 
over twice by computing the mean over adjacent good pixels. 
The pixel sizes are the same in all plots except where noted: 
0.2 dex for both the mass quantites (M, and Mh), 0.3 dex for 
log(l + 85) and 0.15 dex for \og(d prog /Rvir)- Our pixel sizes 
and smoothing scheme were chosen to produce acceptably 
smooth contour lines in order to bring out the main trends. 
Secondary features should not be taken too seriously. 



4.1 M, and Mh for centrals 

We begin with a discussion of the relation between M* 
and Mh. For central galaxies, it has been established with 



high confidence (e.g ., Cacciato et al.ll2009l; More et al. 12009 : 



iMandelbaum et al.lr2o"o9l : iLeauthaud et alll201ll . I2012T ) that 
the relation with the halo mass is somewhat steeper at low 
mass than at high mass. This is a manifestation of the fact 
that the halo mass function and the galaxy mass function 
are considerably different at the low and high mass ends. 
The result is that M* should correlate closely with Mh up 
to around M cr j t ; above this M« begi ns to level off while Mh 
increases (e.g.. lCattaneo et alj|2007l . their Figure 2). So M* 
is a measure of parent Mh for M* below a few xlO lo M0, 
and above that M* serves as a lower limit indicator of Mh. 
In Fig. [4]a. we sho w the relation betw een M, and Mh for 
SDSS central galaxies. lYang et alj (J2008T ) find that galaxies 
above M h ~ 10 12 ' 5 follow a very weak scaling of M* oc M° 22 
while below this mass, M* scales more steeply with Mh, scal- 
ing as M« oc Mj 8 3 for galaxies well below M h ~ 10 125 M© 
(|Yang et al.l 120081 ). We reproduce this result in Fig. |4ji 
(dashed curve) by fitting their Eq. 7 and obtaining very 
similar parameters. This flattening of the M*-Mh relation 
at high masses reflects the fact that massive centrals occupy 
groups with a large range of group membership and group 
masses (since Mh corresponds to group M„), and presents 
an excellent opportunity to separate the effects of quenching 
between M* and Afh as we will see in iJS] 



4.2 Sn and mass for centrals 

As described in ill. II environment density Sn for centrals 
has two modes: the single halo mode measuring field den- 
sity and the cross-halo mode measuring halo mass. Figs. 2J> 
and c show the distribution of SDSS centrals in the <5s-Mh 
and <55-M„ planes. The thick white curve in the top right 
corner of panels b and c marks where the median number 
of group members is six. (More et al., in preparation, give a 
careful treatment of why lines of constant JV ml!m bcrs should 
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Figure 3. The residuals of the SFR division (solid line in Fig. 12b ) 
plotted against the residuals of the rest-frame U — B colour divi- 
sion (solid line in Fig. [T) for the SDSS (0 < z < 0.2). The colour 
bar on the right shows the number density in a pixel smoothed 
over adjacent pixels. The size of a pixel is indicated in the bottom 
left. The galaxies with intermediate SFR lie roughly between the 
two horizontal dashed lines, with the central horizontal solid line 
marking the SFR division. The vertical solid line marks the colour 
division, while the dashed lines roughly delineate the green valley. 
The SFR bimodality mostly corresponds to the colour bimodal- 
ity but 30% of red sequence galaxies lie on the SF sequence while 
only 4% of blue cloud galaxies lie below the SF sequence. Also, 
being in the green valley is not synonymous with transitioning 
from star formation to quiescence. 



run more or less parallel to the thick white lines in Figs. 
|4}) and [5}).) Above this curve, the centrals are mostly in the 
single-halo mode of 85 while below the line, £5 measures dis- 
tances to the next halo. The main distribution of centrals 
in Fig. [4f> (marked by the thin white contours) shows very 
little correlation between 85 and Mh except in the upper 
right corner where single-halo mode centrals are character- 
ized by the highest 8$ and Mh. This region of high Mh, high 
#5 corresponds nicely to the increase of Mh with 85 in the 
upper right corner of Fig. \$p (coloured shading) . We also see 
a deviation from the strong M*-Mh relation in these region 
of high 85 . This simply reflects the fact that higher densities 
in a halo roughly means more satellites, and the summed 
mass of the satellites will begin to matter, compared to the 
mass of the central, in the determination of the halo mass. 



4.3 M*, -Mh and Sn for satellites 

Fig. [SJi shows the M,-Afh relation for satellites. Unlike for 
centrals, there is very little correlation between M, and Mh- 
For small haloes, M* correlates somewhat with Mh in that 
low-mass haloes do not host massive satellites. But the least 
massive satellites occupy host haloes of all masses. 

Environment density is more complex for satellites than 
for centrals. More than half (58%) of SDSS satellites reside 
in groups of six or more members. Fig.© and c show the dis- 
tribution of satellites in the <5s-Mh and <5s-M» planes, again 
with a thick white line marking where the median number 



of group members is six. Above the line, the satellites reside 
in groups of more than six members (the single-halo mode) . 

In the cross-halo mode of <5jv, satellites mostly reside 
in small haloes (M h < 1O 13 M - see Fig. ©) with only a 
few other galaxies. In the cross-halo regions of Fig. © and c 
(under the thick white line), Mh roughly correlates with M* 
(see the coloured shading in these panels). There are at least 
two reasons why M* for satellites is expected to correlate 
with Mh in this regime. 1) We expect a massive satellite to 
live in a massive halo, or else it would be the central. Less 
massive satellites on the other hand are more likely to be 
merging into a low mass halo, since there are many more of 
these than massive haloes. 2) The timescale for dynamical 
friction is inversely pro portional to the ratio of th e satellite 
mass to the halo mass (|Binnev fc Tremainej|l987l ). So very 
large satellites relative to the halo will have disappeared into 
the middle of the halo in a very short time. Therefore, the 
largest remaining satellite must be some fraction of the halo 
mass, but smaller than the ratio of the central mass to the 
halo mass. 

In the single-halo mode, Sn correlates with Mh for satel- 
lites (see the white density contours of Fig. © and the 
coloured shading of Fig. © above the thick white line) just 
as for centrals again because of the halo profile. Since the 
halo profile falls off more slowly in a more massive halo, a 
satellite in a more massive halo will most likely be found near 
the centre of its halo than if it were in a smaller halo, and 
its A?th nearest neighbour would more likely be nearby. For 
a simpler way to understand the Sn-M^ relation for satel- 
lites in the single halo mode, assume a uniform distribution 
of satellites throughout the halo (which is wrong) . Recalling 
that the number of satellites in a h alo is AT sa t oc M£ where a 
is close to 1 (see I Yang et al.ll2008r ). and recalling that in the 
spherical collapse model a virialised halo follows Mh oc R 3 , 
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Therefore en- 



vironment density is a measure of the halo mass following 
something like log(l + <5s) oc | log Mi, which is close to the 
scaling of the SDSS satellites (Fig.©). 

In the single-halo mode for satellites, we also expect 8^ 
to be a measure of radial position within a halo. Fig. ©shows 
85 versus the projected distance d pr oj of SDSS satellites to 
their group centres divided by their group virial radii \R V \ V = 
1.2 x lO 5 (Mh/lO 11 M ) 1/3 pc (|Dekel fc BirnboinJl2009 )]. As 
mentioned in EJ2.1.5I we define D = d pTO j/R v i r . In the single- 
halo mode (above the thick white line which marks where 
the median group membership is six), higher density indeed 
means closer proximity to the group centre. In the cross- 
halo mode, the two quantities do not seem to be related, 
which is what we expect given that 85 measures distances to 
external haloes. (Note that in Fig. ©the pixel size is smaller 
than the usual 0.3 dex for 85 and 0.15 dex for log D in order 
to highlight the valley in the contours that separates the 
single-halo mode from the cross-halo mode.) 

A further complication to this already complex nature 
of 8n is that whether or not a galaxy is in the single-halo 
or cross-halo mode of <5jv depends not only on the number 
of members in its group but also on the limiting magnitude 
of the survey, which is ^-dependent. A group of the same 
mass at z = 0.005, for example, which appears to be in the 
single-halo mode will appear to be in the cross-halo mode 
at z = 0.2 because some of its members will fall beyond the 
magnitude limit. In other words, within the SDSS redshift 
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Figure 4. (a). Stellar mass M, versus halo mass Mj, for SDSS central galaxies (0 < z < 0.2). The dashed curve is a fit to Eq. 7 of Yang 
et al. (2008). (b) Mean stellar mass M* as a function of density log(l + 85) and halo mass M h . (c). Mean halo mass Af h as a function 
of density log(l + 5$) and stellar mass Af*. The thin white contours in each panel mark the weighted number density of points and are 
separated by 0.25 dcx in density. The thick white lines in (b) and (c) mark where the median membership of a galaxy's group is six. 
Membership is greater above this line. The black contours in (b) and (c) follow the coloured shading and are separated by 0.1 dex. The 
number densities in (a) and the means in (b) and (c) were computed in pixels of the displayed size and smoothed over adjacent pixels. 
These panels show that: (1) M* correlates with M^ for low masses, but massive centrals occupy a large range of halo masses; (2) the 
single-halo mode is characterised by an increased 85 and M^ while there is no correlation between £5 and Mj, in the cross-halo mode; 
(3) for low M„, A'/j, increases with M„, and for high M„, M^ also increases with 5s. 
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Figure 5. (a). Stellar mass M* versus halo mass Mh for SDSS satellites (0 < z < 0.2). (b) Mean stellar mass M* as a function of density 
log(l + 85) and halo mass M^. (c) . Mean halo mass Mj, as a function of density log(l + <5s ) and stellar mass M* . The thin white contours 
in each panel mark the weighted number density of points and are separated by 0.25 dex in density. The thick white line in (b) and (c) 
marks where the median membership of a galaxy's group is six. Membership is greater above this line. The black contours in (b) and (c) 
follow the coloured shading and are separated by 0.1 dex. The number densities in (a) and the means in (b) and (c) were computed in 
pixels of the displayed size and smoothed over adjacent pixels. These panels show that: (1) there is very little correlation between M, 
and M h for satellites overall; (2) however M, correlates with M h in the cross-halo mode; and (3) in the single-halo mode <5s correlates 
with Afh as expected in a virialiscd halo with a halo-occupation function that scales with Ar- 



range of 0.005 < z < 0.2, the typical distance to the 5th 
nearest neighbour ranges from 0.34 h" 1 to 4.0 h^ 1 comoving 
Mpc (for 0.005 < z < 0.05 and 0.15 < z < 0.2 respectively). 
In contrast, D does not behave differently in two modes that 
depend on the number of fellow group members. (Recall that 
R v ii depends on Mh which has been corrected for missing 
members). 

Additionally, although both Sn and D suffer from pro- 
jection effects, w e expect the s e effe cts to be worse for 5n 
than for D. The lYang et al.l (|2007l ) group finder includes 
or rejects group members based on a z-window that is de- 
pendent on halo properties (in an iterative way), but 5n 
is computed for all galaxies within a fixed z-window that is 
comparable only to the velocity dispersions of the most mas- 
sive clusters. Since the majority of galaxies are not members 
of such massive clusters, and since the m ajority of satellite s 
in the group catalog are not interlopers (|Yang et alj|2007» . 
we expect 5n to suffer from projection more than D. 



Given these difficulties in interpretating 5m for satel- 
lites, we prefer to use D as a measure of local "environment" . 
But we have included the popular Sn variable in our study 
of quenching in Sj5] in order to compare with previous work. 

For completeness, we replicate Fig. [5f> and c in Fig. [TJi 
and b, except that we put D on the y-axis instead of log(l + 
65). In panel a, we see that the dual-modes of 5 that we 
saw in Fig. [5j> disappear as expected. For example, the two 
regions of high M* in Fig. [5}) become one region in Fig. 0i. 

In Fig. [7J>, the dynamical range of Mh diminishes com- 
pared to Fig. [5J;. This is because low-mass haloes will never 
host very many satellites (nor very massive satellites) and 
#5 for these satellites will almost always be in the cross-halo 
mode; 90% of satellites in haloes of masses 10 12_12 ' 7 Mq live 
in haloes with fewer than six members, so their 85 is forced to 
be low. This is why the lower left corner of Fig.[SJi is coloured 
deep green (low M h ). But these same satellites have no such 
restriction on D. Some of them could reside very near their 



12 Woo et al. 



SDSS: Satellites 
Single — halo mode 




Figure 6. The environment density of SDSS satellites (0 < z < 
0.2) versus their projected distance to their group centres relative 
to their virial radii (d pT0 j / R vlr ) . The thick white line marks the 
dividing line between groups with six or more members (above 
the line) and fewer than six members (below the line). Single- 
halo mode is above this line, cross-halo mode is below. Note how 
the valley in the contour lines tends to follow this line. In the 
single-halo mode, proximity to the group centre is clearly asso- 
ciated with an increase in (5s, while in cross- halo mode, there is 
no relationship. This is expected based on the statistics of these 
relationships. 



group centres even if their 85 is low. Furthermore, if their 
masses are low, they could reside in haloes of all masses. For 
these two reasons, there are no regions in Fig. [7f> with very 
low average Mh. 



5 MAIN RESULTS: QUENCHING AND HALO 
ENVIRONMENT 

In the previous section we explored the relations between 
stellar mass, halo mass and environment density for satellite 
and central galaxies in the SDSS. We showed that for central 
galaxies, M* correlates with Mh for low masses, but massive 
centrals of a given stellar mass occupy a large variety of halo 
masses. We also showed that the single-halo mode of 85 for 
centrals is characterised by increasing Mh with increasing 
85 while in the cross-halo mode, there is little correlation 
between 85 and Mh. For satellite galaxies, there is a weak 
correlation between M, and Mh in the cross-halo mode only. 
In the single-halo mode, 85 for satellites correlates with 
Mh as expected for virialised haloes with a halo-occupation 
function that scales with Mh. Since we now know how these 
quantities relate to each other, we are in a position to explore 
quenching as a function of these quantities. 



5.1 Quenching for centrals 

Fig. [8] repeats the panels of Fig. [4] (M* and Mh, <5s and M h , 
and 85 and M,) but now shows the fraction of quenched 
galaxies as the third parameter represented by the colour 
gradient. Fig. [SJj is the key result of this paper. It shows 
that the quenched fraction appears to increase more strongly 



with Mh at fixed M* than with M* at fixed Mh for central 
galaxies. The contour lines of the quenching to first order 
follow Mh, not M*, i.e., they are roughly vertical above 
M* ~ 10 10 ' 4 Mq. The quenched fraction increases by less 
than ~ 0.1 over 0.8 dex, or 40% of the range in M* (at con- 
stant Mh), but increases by almost 0.35 over 40% of the Mh 
range (at constant M*). This result seems to prefer halo- 
quenching models over mechanisms that are more closely 
tied to M* such as the QSO mode of AGN feedback. Note 
that the halo environment heated by a virial shock is more 
likely to host the radio mode of AGN (jDekel fe Birnboiml 
2006) because the coupling of the AGN energy to the halo 
gas is more efficient when the gas is hot. 

To second order, Fig. [8)2 shows a small dependence of 
quenching on M* at constant Mh for massive haloes. This 
can be explained in the following way. After a virial shock 
forms in a massive halo, residual cold gas may still form 
stars inside the galaxy after the shock has been triggered. 
An M* dependence could then naturally arise if at fixed 
Mh, a central with higher stellar mass had converted more 
of any residual gas into stars, exhausting its gas supply, and 
making it more susceptible to quenching. However, this M*- 
dependence of the quenching seems to be a secondary effect 
compared to the dominant trend of quenching with Mh as 
seen in Fig. [HJx. 

Fig. [8j> shows quenching in the S^-Mh plane. It re- 
veals that the quenched fraction increases with Mh and only 
weakly with 85 . At constant Mh, the quenched fraction in- 
creases at the most by 0.2, and for most masses only by 0.1, 
over two orders of magnitude in 85. However at constant 55, 
the quenched fraction increases by as much as 0.4 over two 
orders of magnitude in Mh . 

When we examine quenching in the 8$-M* plane 
(Fig. [HJ;), especially at large 85, we see that the quenched 
fraction for centrals depends on both 85 and M*. The in- 
crease of the quenched fraction with M* is consistent with 
the increase of M h with M„ below M* ~ 10 n M Q (see 
Fig. |4};). Above this mass, the quenched fraction increases 
with £5 in a way similar to the increase of Mh with 85 in 
Fig. \4fc. Thus, apart from the low M* and high 85 region, 
the quenching pattern in Fig. [5J; is explainable by a basic 
trend of quenched fraction with Mh. In Sj?]we discussed that 
the Mh pattern in Fig. 03 is the result of the interplay be- 
tween the two modes of Sn and the weakening of the M*-Mh 
relation for massive centrals. Thus the quenching pattern in 
the ($5 — M* plane is con sistent with a corr elation of quench- 
ing with M h . Note that IPeng et all (|201lr i did not find the 
85 dependence of the quenched fraction at high M* mostly 
because they denned quenching by colour and not by SFR - 
more on this in £15.31 

The dependence of quenching on 85 at low to interme- 
diate Mh and at high 85 (see Fig. [8j), also reflected in c at 
low M t ) is not a neccessary trend in the simple theoretical 
framework of halo quenching that we have been advancing. 
This may possibly be attribu ted to misclassifie d satellites 
(an explanation suggested by IPeng et al.l I2011T ) . However, 
we can offer another explanation that takes the classification 
at face value and appeals to the halo quenching mechanism. 
Halo masses on the order of Merit may be capable of sustain- 
ing a stable virial shock without actually containing one. It 
has been seen in cosmological simulations that, once the halo 
mass is ~ M cr it, the trigger for the formation of a shock that 
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Figure 7. (a) Mean stellar mass M* as a function of relative group distance D = d pro j/Rvir and halo mass M^. (b). Mean halo mass M^ 
as a function of relative group distance D and stellar mass M, . The thin white contours in each panel mark the weighted number density 
of points and are separated by 0.25 dex in density. The black contours in follow the coloured shading and are separated by 0.1 dex. The 
means were computed in pixels of the displayed size and smoothed over adjacent pixels. These panels show that: (1) the D-M^ plane 
does not exhibit dual-mode behaviour as in Fig. [5}>; (2) the dynamical range of mean M^ in the D-M t plane is diminished compared to 
the 5&-M* plane. 



quickly expands to the virial radius is typically an occasional 
minor merger. Small centrals that have a high 8n must have 
one or more external haloes nearby, and are therefore more 
likely to capture loose satellites of these neighbouring haloes, 
which may trigger the virial shock that causes quenching. 
This explanation predicts a higher quenched fraction at the 
highest values of 5$ in halo somewhat below M CT n , as is ob- 
served. 

An astute reader may notice that the quenched frac- 
tion seems to reach lower levels when using M* (Fig. [§);) 
than when using M h (Fig. [Hf)). These lower levels of quench- 
ing at low M* compared to low Mh reflect the fact that 
the quenching contours in Fig. [8ji are close to horizontal 
for M. < 1O 1O ' 4 M and M h < 1O 12 ' 5 M . This is the regime 
where Mh-dependent quenching is not expected to be im- 
portant (since Mh is below or around M cr it)- It is here that 
residual M*-dependent quenching becomes important, and 
this regime will be explored further in ij6.ll 

A correlation between the quenched fraction and halo 
mass at fixed stellar mass, for centrals, is expected on the- 
oretical grounds bas ed on the critical mass for viria l shoc k 
hea ting, as inlBirnboim fc Dekell J2(X)3) , iKeres et all (|2005h , 
and lDekel fc Birnboiml (J200fj). We note, however, that a de- 
tected correlation does not necessarily imply causation, and 
that one should apply caution because the halo mass and 
the stellar mass are not independent of each other. Indeed, 
at low masses, A/* < 10 n M Q and M h < 1O 12 M , halo mass 
and stellar mass are expected to grow in concert, while at 
larger masses the stellar ma ss tends to a constant as the 
halo mass keeps g rowing fe.g.. lCoorav fc Milosavlievial2005l ; 
lYang et al.l 120121 ). Can this introduce an apparent correla- 
tion of the quenched fraction with Mh at a fixed M, in 
an alternative scenario where the "true" cause of quench- 
ing is associated with M* and not with M^l It can clearly 
introduce a general trend, where the average halo mass of 
quenched galaxies is larger than for un-quenched galaxies, 



but it is not expected to introduce a continuous gradient 
of the quenched fraction with Mh at a given M* all the 
way to 1O 15 M0, and no correlation of the quenched fraction 
with M* at a given Mh, as in Fig. [H}a. We conclude from 
Fig. [HJz that the direct contribution of halo mass to quench- 
ing seems to be more fundamental than the contribution of 
stellar mass. 



5.2 Quenching for satellites 

Fig. [9] shows the equivalent of Fig. [8] but for satellites. In 
all three panels of Fig. [9] the quenched fraction for satellites 
increases with all three quantities M^ , M* and <5g suggesting 
that these three projections are not the most helpful way to 
understanding the main drivers of quenching for satellites. 
As argued in Sj4] Sn for satellites is a complicated quantity, 
whose properties change according to the choice of TV, on 
limiting magnitudes. Furthermore, M* for satellites is ex- 
pected to correlate with their position in their host halo, 
e.g., because dynamical friction causes massive satellites to 
spiral inward faster than less massive ones. Moreover, less 
massive satellites suffer from colour-dependent incomplete- 
ness. A mass-limited cut on the satellite sample which would 
m ake it complete for red galaxies (following the cut applied 
bv lDutton et al.ll201ll : M, > 10.7+log(z/0.1)) includes very 
few galaxies smaller than M* ~ 1O 1O,3 M and excludes more 
than half the satellite sample. 

In light of these difficulties in interpreting Fig. [5] we 
complement Figs. [9j> and c with Figs. 110b and b, replacing 
($5 on the y-Sixis with the simpler projected radial distance 
from the centre (D = dp ro j/-Rvir). (Recall from Fig. [6] that 
Sn correlates with D in the single-halo mode of 5m-) 

Fig. 110b shows that the quenched fraction for satellites 
increases with both halo mass and with proximity to the 
group centre. The quenching trend with radial distance is 
strongest for distances of greater than ~0.1-0.2i? v i r . Below 



14 Woo et al. 



SDSS: Centrals 




SDSS: Satellites 



12.0 12.5 13.0 13.5 14.0 14.5 15.0 
log M h /M 



o 




12.0 12.5 13.0 13.5 14.0 14.5 15.0 
log M h /M 



□ 




n 



n 



(c) 



9.0 9.5 10.0 10.5 1 1.0 1 1.5 12.0 

og M„/M 

Figure 8. The mean quenched fraction as a function of M, and 
M h (a), of ^5 and M h (b), and of £5 and M* (c) for central galaxies 
in the SDSS (0 < z < 0.2). The thin white contours in each panel 
mark the weighted number density of points and are separated by 
0.25 dcx in density. The thick white line marks where the median 
membership of a galaxy's group is six. Membership is greater 
above this line. The black contours follow the coloured shading 
and are separated by 0.05. The means were computed in pixels 
of the displayed size and smoothed over adjacent pixels. To first 
approximation, the quenched fraction for centrals correlates with 
Mi! more strongly than with M* or 5s as shown by the vertical 
black contours in panel a and b. Quenching in c roughly follows 
mean M h (see Fig.[4j:). 
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Figure 9. The mean quenched fraction as a function of M* and 
M h (a), of #5 and M h (b), and of £5 and M„ (c) for satellite galax- 
ies in the SDSS (0 < z < 0.2). The thin white contours mark the 
weighted number density of points and are separated by 0.25 dcx 
in density. The thick white line marks where the median mem- 
bership of a galaxy's group is six. Membership is greater above 
the thick line. The black contours follow the coloured shading 
and are separated by 0.05. The means were computed in pixels of 
the displayed size and smoothed over adjacent pixels. Quenching 
for satellites correlates with all three quantities M h , M t and 5$. 
None of these panels contains regions of very high mean quenched 
fraction, indicating that, in contrast to centrals, the combination 
of M„, M h and 5$ is not optimum for expressing the quenched 
fraction of satellites. 
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Figure 10. The mean quenched fraction of SDSS satellites (0 < z < 0.2) as a function of the projected distance to the group centre 
relative to the group virial radius (d pro j / R v ir) and the halo mass Mh (a) and stellar mass (b). The thin white contours in each panel 
mark the weighted number density of points and are separated by 0.25 dex in density. The black contours follow the coloured shading 
and are separated by 0.05. The means were computed in pixels of the displayed size and smoothed over adjacent pixels. Quenching for 
satellites increases with M h and decreases with D = d pro j/R v i I . Quenching also increases with M t at large D where M^-qucnching is 
weakest and where the satellites likely fell in recently, suggesting a quenching dependence on the sub-halo. Satellites closer to the group 
centre than ~ 0.2iJ v i r likely fell in less recently, and their quenched fraction depends strongly on host Mh and almost not at all on 
satellite M*. 



this, the motions of satellites and/or projection effects may 
smear out the quenching trend with D, causing quenching to 
appear to increase only with Mh. However we do not expect 
projection effects to be strong, especially in massive haloes 
which contain many satellites at small radial distance so the 
signal-to-noise for small D is higher. 

Fig. 110b shows that the quenched fraction also increases 
with satellite stellar mass M, for constant D. The M» de- 
pendence is strongest for satellites at the farthest distances 
from the centre, where the Mh dependence is weakest (but 
even here, the Mh-dependence of the quenched fraction is 
not that weak, increasing by more than 0.3 over 2.5 dex). If 
we assume that these galaxies at large D have only just fallen 
into a larger halo, their quenching dependence on M* is con- 
sistent with a quenching dependence on the mass of their 
own halo before infall, now sub-halo, since Mh-dependent 
quenching operated on them while they were still centrals 
(Fig. [8jz). The M, dependence in Fig. 110b virtually disap- 
pears for satellites closer to the centre than ~ 0.2i? v i r , where 
the Mh dependence is the strongest. If these satellites have 
resided in their host haloes for some time, the quenching 
processes which operate in the host halo and which depend 
on the host Mh will eventually affect the satellites so that 
their quenched fraction will depend strongly on host Mh 
rather than on the satellite M* as we see in Fig. 1101 

Fig. 1 10b succeeds in identifying the properties of satel- 
lites that produce high quenched fractions of > 0.95, namely 
high Mh and close proximity to the halo centre. None of the 
panels in Fig. [9] or Fig. 110b produce regions of the plot with 
such high quenched fraction. This seems to suggest that it is 
a combination of halo mass and distance to the halo centre 
that governs quenching in satellites. 

In terms of dynamical range of quenching, however, the 
combination of D and Mh produces a range as broad as the 
pair <5 5 and M h (Fig. [9};) or even the pair 85 and M* (Fig. [9}:). 
However, quenching trends with 85 are difficult to interpret 



due to the dual nature of 8. Which mode of 8 a satellite 
belongs to depends on redshift as discussed earlier, in that 
a galaxy in the single-halo mode at low z will appear to be 
in the cross-halo mode at high z because some members of 
its group will drop out of the magnitude limit. Furthermore, 
these satellites with spuriously low 85 will likely be bluer on 
average than the true population, especially at lower mass, 
since low mass red galaxies are selected against at higher 
z. Indeed, when we checked the sub-sample with z < 0.05, 
satellites with -1 < 8 5 < and 9 < logM*/M Q < 10.5 had 
quenched fractions 0.1 higher than the whole redshift sample 
in these 85 and M* ranges, whereas the quenched fraction in 
the D-Mh plane did not change. For these reasons, quench- 
ing trends with 85 may be somewhat exaggerated, and D is 
the preferred predictor of quenching. 

As for whether Mh or M* is better at predicting quench- 
ing for satellites, Fig. |5}j seems to indicate that they are 
equally good. However, quenching seems to follow average 
Mh better than average M* . Specifically, the quenching con- 
tours in Fig. [SJ: roughly follow the Mh contours in Fig. [SJ;, 
while the quenching contours in Fig. |9}> do not follow the 
M* contours of Fig. [SJ>. Similarly, the quenching contours 
in Fig. 110b do not follow the M, contours in Fig. [TJz (ex- 
cept at high D) while the quenching contours in Fig. 110b do 
somewhat follow the Mh contours in Fig. [7J) (except at high 
D). (Quenching at the high D as we argued above may be 
following sub-halo mass more than host halo mass.) 

In summary, we have found that the combination of 
halo mass and distance to the halo centre is an excellent pre- 
dictor of quenching for satellites, producing high quenched 
fractions of > 0.95. Combinations involving 55 and M* also 
produce large dynamic range of quenching, but this may be 
due to a lucky combination of redshift effects and sub-halo 
effects. 

The trend of the quenched fraction with radial distance 
from the group centre is consistent with external media- 
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nisms of quenching such as the halo quenching model, where 
virial shock heating is most efficient for large halo masses 
and in the centres of haloes. In such a scenario, quench- 
ing is expected to be a function of the density of the hot gas 
which increases with decreasing radial distance to the centre 
as well as with Mh. This trend of quenching with proximity 
to the centre is also consistent with other processes such as 
strangulation, tidal stripping and ram pressure stripping. 

5.3 Comparisons with previous work 

Our finding that quenching increases primarily with Mh 
rather than M*, and is also correl ated with D for s a tellites , 
is consistent with the findings of IWeinmann et al.l (|2006al ) 
who find that halo mass is more important than luminosity 
in predicting galaxy properties (and improves on their result 
as M* is a more intrinsic internal property than L). 

Our study improves on that of lKimm et al.l (|2009f ) who 
show similar plots as Fig. [8ji and Fig. [9ji and could not de- 
cide which of M* and Mh was more important for quench- 
ing in centrals. For satellites they argued that M* and Mh 
were equally important for quenching. We have improved on 
their analysis by combining M»- and ^-dependent quench- 
ing into a dependence on D and Mh- For centrals, we have 
improved on their result by smoothing the quenched frac- 
tion over the M,-Mh plane in order to discern the overall 
trend of quenching with Mh- The slight M» -dependence of 
the quenched fraction for c e ntrals with high Mh is similar to 
that found in iKimm et al.l i|2009l ). but we have shown that 
the the overall quen ching trend with M h is stronger. 

Note also that IKimm et al.l (|200gf ) use stellar masses 
that were computed using the models of iBell et al.l (|2003l ) 
which can underestimate M* for dusty, star-forming galaxies 
(for a typical dust model with attenuations of 1.6 and 1.3 in 
the a and r bands, M» will be underestimated by 0.2 dex: 
see IBell et al.l 120031 ). Underestimating M* for star forming 
galaxies but not for quiescent galaxies has the effect of exag- 
gerating quenching trends with M, at fixed Mh- In contrast, 
the values of M» that we use are estimated by SED fitting 
taking dust into account (see the references in H2.1.4|l . 

One may w orry that the gr oup stellar mass es of 
I Yang et al.l (|2007l ). computed using IBell et all (120031 ). also 
result in exaggerated quenching trends with M h - However, 
we do not expect this effect to be large especially above 
Mh ~ 1O 13 M0 where the centrals are more than 70% 
quenched and their masses contribute less than 75% of their 
group mass. To test the effect of dust on the Mh estimates, 
we performed a rough rescaling of the Yang et al. Mh values 
to reflect group masses based on Brichmann et al. M« esti- 
mates, and find that the quenched fraction increases more 
than twice as much with Mh at fixed M„ than with M* at 
fixed Mh over the same range of mass. This is despite the 
uncertainties in Mh being much larger (> 0.3 dex) than the 
typical uncertainties in M„ (< 0.05 dex). So, after consider- 
ing dust, our result that Mh predicts quenching more than 
Mt, appears to be robust. (More on dust below.) 

Our result that quenching of satellites correlates with 
Mh an d anticorrelates with D i s co nsistent with t he re - 
sults of iBlanton fc Berlindl J2007I ) and lHansen et all l|2009l ) 
(who use group luminosity and group richness rather than 
Mh). The decreased quenching with D also agrees with 
several other studies of SFR or the blue/red fractions 



within gro ups and clusters (iGomcz et al 
2004J; iTanaka et al.l liooi; feines et all 



20031 ; iBalogh et al 



2005; 



Haines et al 



20071; I Wolf et al.ll2009l:lvon der Linden et al.l2010h. (See also 



Diaferio et al. 



] l200ll ; IWeinmann et al.l 
20101 for semi-analytic model recipes that reproduce these 



trends). 

Ou r results also agree with those of IWetzel et al] 
|2012al ) who find that qu enching for satellites c orrelates with 
M h and D at fixed M„. IWetzel et all (l2012bl ) also demon- 
strated that the distribution of specific SFR for satellites 
can be ex plained by cons i dering the satellites' time since 
first infall. |Pe Lucia et al] (|2012l ) also showed using a semi- 
analytic cosmological model that satellite distance from the 
group/cluster centre is correlated with the time of accre- 
tion. This is consistent with our finding that the satellite 
quenched fraction increases with M, at higher D values. 
Since these satellites recently fell in, their quenched fraction 
reflects their sub-halo mass rather than their h ost halo mass. 
But so me time after the first infall (2-4 Gyr: IWetzel et al] 
l2012al ; 5-7 Gvr: lDe Lucia et aljliorl ), halo quenching takes 
over so that quenching for satellites closer than ~ 0.2i? v j r 
depends st rongly on their host h alo mass. 

While lBamford et al.l |2009t ) found that the red fraction 
of galaxies increases with M* and environment (as measured 
by Sn and the distance to the nearest galaxy group), they 
also found that the red fraction does not vary strongly with 
group mass (see their Fig. 13). This at first sight seems to 
contradict our result that quenching depends on Mh (de- 
rived from group mas s). However, these authors use the 
C4 Cluster Catalog of iMiller et~atl (|2005l ) which does not 
include isolated galaxies. The membership of each group 
ra nges from 10 to more than 200 galaxies so that Fig. 13 
of iBamford et al.l (|2009l ) is dominated by satellites. As we 
showed in Fig. 110b . both Mh and D are needed to predict 
quenching for satellites. Most satellites lie at 0.3-li? v ir from 
their group centres (white contours of Fig. 110b ) which is 
where the Mh-dependence of quenching is weakest. There- 
fore, measuring the quenched fraction of satellites as a func- 
tion of Mh , regardless of D, will result in a weak Mh depen- 
denc e of quenching. This also se ems to be the main reason 
why Ivan den Bosch et al] (|2008l ) find only weak correlation 
between satellite colour and Mh and between satellite colour 
and D. 

Our analysis imp roves on the interpretation of 
IPeng et al.l (|2010l . l201ll ) who find that mass determines 
quenching for both centrals and satellites, while 85 also gov- 
erns quenching for satellites. They do not distinguish be- 
tween stellar mass and halo mass for centrals (for satellites, 
see below). While quenching for centrals apparently depends 
on M, (Fig. [8)3), we have found that quenching correlates 
more strongly with Mh at fixed M* than with M* at fixed 
Mh (Fig. [8J3.) . We also related the £5 and M* dependence 
of quenching for satellites to a quenching relation with D 
and Mh and argued that these are the preferred predictors 
of quenched fraction. 

However. IPeng et al] |201ll ) find that for satellites, M,- 
quenching is more important than Mh-quenching at a fixed 
65, namely for 1.0 < log(l + S$) < 1.3. But this is not the 
whole quenching picture. Using Fig. [6] this narrow range 
of 5s corresponds to a log D of about -0.4 which is in the 
region of Fig. [10] where satellites only recently fell in. Here, 
quenching relates to the sub-halo and the host halo mass 
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Figure 11. The mean red fraction as a function of M* and Mh for centrals (a) and satellites (b) in the SDSS (0 < z < 0.2). The thin 
white contours in each panel mark the weighted number density of points and are separated by 0.25 dex in density. The black contours 
follow the coloured shading and are separated by 0.05. The means were computed in pixels of the displayed size and smoothed over 
adjacent pixels. To first approximation, the red fraction for centrals correlates with M* more strongly than with Mj, for both centrals 
and satellites. The red fraction behaves differently from the quenched fraction (Fig. [HJi, |5Ji) because red galaxies which are still forming 
stars are on average more massive than the star-forming population. 



will not matter much. But we demonstrated that quenching 
of satellites strongly increases with Mh (and almost not at 
all with M*) for satellites closer to the group centre than 
~ 0.2.R v i r where we expect the time since infall to be long. 

5.4 SFR vs. Colour 

Lastly, we caution that measuring quenching by colour 
rather than by SFR causes M*-dependent quenching to 
appear stronger and Mh-dependent quenching to appear 
weaker than they actually are. The reason is that colour 
as a quenching proxy is severely contaminated by dust. As 
described in SJ31 one third of the red sequence galaxies are in 
fact not quenched but are dusty star-formers. These predom- 
inantly lie on the massive end of the SF sequence (but share 
a similar mass distribution to that of the red sequence). 
Thus, dividing the galaxies by colour removes many of the 
massive galaxies from the "star forming" sample, leaving an 
overall less massive "blue" sample (by about 0.4 dex). (This 
same exercise does not render the "red" sample with a signif- 
icantly different mass distribution as the "passive" sample). 
Therefore, the red fraction will appear to correlate with M* 
more than the quenched fraction. Furthermore, since M* and 
Mh do not correlate for satellites, and only very weakly for 
centrals at the massive end, the red fraction will appear to 
corr elate more w i th M* than with Mh. And this is the reason 
why IPeng et al.l (|201lT ) missed the significant (^-dependent 
quenching for centrals at the massive end (Fig. |8p) that we 
showed follows the increase of Mh with 85 at these masses 
(Fig. Sb) , reflecting the stronger trend of quenching with Mh 
than with M*. This is also the reason why they are led to 
find that Mh is unimportant to the red fraction for satellites. 
To illustrate this effect of dust, we show the red frac- 
tion for centrals as a function of M» and Mh in Fig. Illb . 
reproducing Fig. [SJi, but using the red fraction as the third 
parameter denoted by the coloured shading. The red frac- 
tion increases almost entirely in the M* direction in contrast 
to the quenched fraction (Fig. [SJz). Fig. Illb is the analogous 



figure for satellites. Comparing Fig. Illb with Fig.[9jx, the red 
fraction increases more strongly with M* than with Mh, but 
the quenched fraction seems to increase equally str ongly for 
both masses (which agrees with lKimm et al.ll2009T ). In both 
panels of Fig. 1111 for constant Mh, the galaxies with higher 
M* will have a higher fraction of dusty star-formers than 
those with lower M», so that correcting for dust will make 
the quenching lines become more vertical. The striking dif- 
ference between the red fraction and the quenched fraction 
demonstrates that dust can mislead one's interpretation of 
the primary drivers of quenching. 

Our result that Mh is more important than M* for 
quenching represent s a significant physical difference from 
the interpretation of lPeng et al.l (|2010l . 120111 ) , but their ana- 
lytic represention of quenching might still be a useful math- 
ematical tool. These authors proposed that the dual depen- 
dence of the quenched fraction on 85 and M* is separable, 
with 85 governing the quenching of satellites and M* gov- 
erning the quenching of centrals. Our results improve their 
model by substituting M* with Mh with no dependence of 
quenching on 85 for centrals (the £5 -dependence at constant 
M» or Mh is very weak except at the highest 85 - see Fig. [8}? 
and c and the discussion in J)5.1|) . For satellites, one may 
substitute 85 and M* with D and Mh assuming that D- 
quenching and Mh-quenching are separable. 



6 THE EVOLUTION OF THE RELATIONS 
BETWEEN STELLAR MASS, HALO MASS 
AND DENSITY 

The relations between stellar mass, halo mass and envi- 
ronment density as described in Sj3] are expected to change 
with redshift. The most important physical reason for such 
a change is that the halo mass fu nction is expected to grow 
between z ~ 1 to z ~ 0. The iPress fc Schechterj (|l974T ) 
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cordance cosmology. The number density of haloes with 
M h = 10 13 ' 6 , 10 14 ' 5 , and 10 15 M Q grows by a factor of 2, 
6, and 30 from z = 1 to z = (using the halo mass func- 



ti on of Tinker et al.l 20081 and the mass transfer function 
of lEisenstein fc Hul |l998J)). Therefore many more haloes at 
z ~ 1 than at z ~ will be below Merit. In the AEGIS 
field, 85% of the galaxies at 0.75 < z < 1 live in haloes less 
than M h = 10 12 ' 5 M Q , in contrast to only 39% of the SDSS 
galaxies below this mass. Thus at z ~ 1 we expect to see 
if any Mh-independent quenching becomes dominant in the 
absence of Mh-dependent quenching. 

A second difference between z ~ 1 and z ~ is due 
to the brighter magnitude limit at the higher redshift re- 
sulting in fewer satellites observed at z ~ 1 than at z ~ 0. 
The weighted satellite fraction for the AEGIS sample in the 
redshift bin 0.75 < z < 1 is 10% compared to 30% in the 
SDSS above M* = 10 1 Mq. Since the number of satellites in 
the AEGIS sample is small (162) any conclusions concerning 
quenching for these satellites would be dubious and so we 
choose not to address satellites for this redshift. 

In order to see how the relations between Mh, M* and 
6n (N = 3 for AEGIS; see §2]) might change at z ~ 1 due 
to reduced halo masses, we performed the same analysis on 
the AEGIS centrals as we performed on the SDSS centrals. 
(Note the following figures depicting a third quantity as a 
function of the x- and y-axes are computed using the same 
smoothing as in the SDSS with m — 2; see i|4]) 

Fig. [12] is the high-z (0.75 < z < 1) analogue to Fig. 0] 
(i.e., the M,-M h , S 3 -M h and <5 3 -M* planes). Fig. [I2ji shows 
that the stellar mass scales with the halo mass for most of the 
Mh range, and starts to flatte n around Mh ~ 10 12 5 Mq. We 
show the best fit of Eq. 7 of I Yang et al.l 1 20081 ) for central 
galaxies in Fig. 112b as the green dashed curve. The high- 
end slope of 0.21 is similar to that of SDSS for centrals 
(blue dotted curve of Fig. I12[) . but the knee point is 0.5 
dex higher, placing the upper part of the M*-Mh relation 
for AEGIS higher than that of the SDSS. This is consistent 
with halo growth since z ~ 1 around centrals that do not 
grow much in M». 

Our M»-Mh relation at z ~ 1, and its normalization rel- 
ative to the Mt-Mh r elation at z ~ 0, are in g ood agreement 
with the findings of iLeauthaud et al.l (J2012J) who measure 
this relation for the COSMOS survey at z ~ 0.9 in compar- 
ison with th e SDSS. Our relatio n at z ~ 1 is also consistent 
with that of iMoster et al.l Q2010J) quantitatively, however its 
relative normalisation compared to the M«-Mh relati on at 
z ~ differs from their findings. IMoster et al.l (|2010P ) find 
that the upper portion of the z ~ 1 M*-Mh relation roughly 
coincides with the relation at z ~ while the lower portion 
shifts to lower M* at z ~ 1 at fixed M h relative to the re- 
lation at z ~ 0. ILeauthaud et al.l (J2012J) speculate that this 
discrepancy is due to underestimated errors on the part of 
Mostcr et al. (2010). Further discussion on the origin of the 
discrepancy is beyond the scope of this paper. 

The number of observed group members is almost al- 
ways below N = 3 for centrals (fewer than 1% of centrals 
reside in groups of more than three observed members, and 
these live in the most massive haloes). Therefore we expect 
to see in AEGIS only a muted form of the two-component 
distribution in 8m vs. Mh that we saw for centrals in the 
SDSS data (Fig. |4jj). Most of the galaxies being in the cross- 
halo mode, S3 should not correlate with Mh except at the 



very highest masses. As expected, Fig. 112b shows that the 
distribution of AEGIS centrals shows no correlation between 
S3 and Mh, except that the centrals in the most massive 
haloes never reside in low densities. For the same reason, 
the coloured shading of Mh in Fig. 112b runs with M, and 
not with S3, except perhaps weakly at the very massive end. 

6.1 Quenching at high-z 

As we have argued, halo quenching via virial shock seems to 
dominate the quenching of galaxies in haloes above Mh ~ 
10 Mq, i.e., those that are large enough to support a virial 
shock. However, haloes smaller than this are less likely to 
support a stable virial shock and such haloes would be ideal 
for testing whether non-halo-dependent quenching mecha- 
nisms operate when virial shock heating is relatively unim- 
portant. As described above, the majority (85%) of galaxies 
at z ~ 1 reside in haloes smaller than Mh ~ 10 12 ' 5 Mq (cor- 
responding to Af* ~ 10 11 Mq for centrals - see Fig. 112b ) 
while the stellar mass function of galaxies has remained rel- 
ative ly constant over time since z ~ 1 (e.g., iBundv et al.l 
|2006J). While the halo quenching model predicts very low 
quenched fractions at all masses below Mh ~ 10 12 ' 5 M©, the 
behaviour of the quenched fraction in this regime will indi- 
cate the strength of other quenching mechanisms when halo 
quenching is expected to be less important. 

Fig. 1131 shows the same panels as Fig. [12] (i.e., central 
galaxies at z ~ 1) with the quenched fraction as the third 
parameter denoted by the coloured shading. Panel a shows 
that the quenched fraction depends only on M, (horizontal 
red lines) for low-mass centrals but begins to depend also on 
Mh at higher masses (tilted red lines). The increase of the 
quenched fraction with M* right below Mh ~ 1O 12 ' 5 M0 re- 
flects a similar behaviour at 2 ~ (see Fig. [HJj at the lowest 
masses) and points to an M*-dependent quenching mecha- 
nism that operates where halo quenching is not expected to 
dominate. 

Note that our AEGIS sample is co mplete only down 
to M, ~ 10 10 ' 7 M Q (|Noeske et all 120071 ) . and the missing 
galaxies below this mass are likely to be red or quiescent 
and of lower mass. Therefore, for fixed Mh in Fig. 113b . re- 
gions of lower M„ in the incomplete region may actually be 
more quenched than observed. In other words, the quenching 
contours in Fig. 113b may in reality be more tilted than hori- 
zontal. However we do not ex pect the effect to be very large 
since there is evidence (e.g., iBell et aLll2004l ; iBrown et al.l 
120071 ; iFaber et al.ll2007l ) that there are fewer red galaxies at 
z ~ 1 to be missed at all. 

In panel c of Fig. 1131 the quenched fraction increases 
with M* for centrals less massive than M» ~ 10 1 Mq, but 
begins to increase also with log(l + S3) above this mass. For 
these galaxies the increase of the quenched fraction with S3 
coincides with an increase in Mh (Fig. 112b ) similar to what 
we saw for the SDSS centrals (JSTTJ. 

In order to ensure that the M*-dependence of quench- 
ing we observe here is not due some artefact of the group 
finding process, i.e., contamination by satellites, we tested 
the AEGIS mock catalogues in the following simple manner. 
For each mock catalogue, we assumed a mass-independent 
quenching such that 20% of all centrals and 80% of all satel- 
lites are quenched regardless of M» or Mh. We then com- 
puted the quenched fraction in 10 bins of M, and Mh for 
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Figure 12. (a). Stellar mass M* versus halo mass Mj, for the AEGIS (0.75 < z < 1). The green dashed curve is a fit to Eq. 7 of Yang 
et al. (2008). The blue dotted line is the analogous fit for the SDSS. (b) Mean stellar mass M« as a function of density £3 and halo mass 
M^ . (c) . Mean halo mass Mj, as a function of density log( 1 + S3) and stellar mass M» . All panels show only central galaxies in the AEGIS 
sample. The thin white contours in each panel mark the weighted number density of points and are separated by 0.1 dex in density. The 
black contours in (b) and (c) follow the coloured shading and are separated by 0.1 dex. The number densities in (a) and the means in 
(b) and (c) were computed in pixels of the displayed size and smoothed over adjacent pixels. The M^-M^ relation for massive centrals 
at z ~ 1 lands higher than that of the centrals at z ~ 0, consistent with halo growth since z ~ 1. £3 and Mj, do not correlate except that 
centrals of massive haloes do not live in sparce environments. 
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Figure 13. The mean quenched fraction as a function of M* and Mh (a), of £3 and Mh (b), and of S3 and M* (c) for central galaxies in 
the AEGIS (0.75 < z < 1). The thin white contours mark the weighted number density of points and are separated by 0.1 dex in density. 
The red (dark) contours follow the coloured shading and are separated by 0.05. The means were computed in pixels of the displayed 
size and smoothed over adjacent pixels. Below Afj, ~ 10 12 ' 5 Mq, quenching for centrals at z ~ 1 increases primarily with M*. Above 
this mass, quenching increases with £3, corresponding with the increase of Mh with £3 in Fig, ll2fc just as for the SDSS centrals. But the 
overall quenched fraction remains low. 



centrals as determined by the group finder. The Pearson cor- 
relation coefficients between the quenched fraction and mass 
were on average 0.02 ±0.08 for the quenched fraction vs. M* 
and 0.14±0.08 for the quenched fraction vs. M h . For the real 
AEGIS data, the Pearson correlation coefficients are about 
0.9 between the quenched fraction and both A/* and Mh- In 
other words, satellite contamination in the quenching trends 
is too low to explain the quenching correlations with mass 
seen in Fig. 1131 

Regardless of these trends of the quenched fraction, 
the total quenched fraction remains low at all masses. The 
weighted quenched fraction for all the centrals at (0.75 < 
z < 1) is 0.11. Since the DEEP2 .R-band magnitude limit 
corresponds to a lower M, for blue galaxies than red galax- 
ies, we compare these weighted quenched fractions with 
those of the SDSS after applying a completeness cut that 
samples the same complete region o f colour-mass space out 
to z = 1 (refer to lGerke et al.ll2007l for details). The colour- 



mass-complete weighted quenched fraction is 0.41 for cen- 
trals, almost four times those of the AEGIS galaxies. This 
dearth of quenched galaxies at z ~ 1 compared to z ~ 
is consistent with the growth of the red sequence over the 
last ~ 7 Gyr (e .g., iBell et all |2004| ; iBrown et all |2007l : 
iFaber et al.ll2007l ) . Since the halo mass function is predicted 
to have grown an order of magnitude in mass since z ~ 1, 
the quadrupling of the colour-mass-complete quenched frac- 
tion from z ~ 1 to z ~ can naturally be explained by halo 
mass dependent quenching. 

Our result that M* -dependent quenching dominates 
over Mh- and 6*3- dependent quenchin g for centrals at z ~ 1 
is consistent with the predictions of iMcGee et al.l (J2009T) . 
These authors used a semi-analytic model to predict that 
the fraction of "environmentally affected" galaxies should 
be very low at z~ 1.5. 

Although Mm -dependent quenching seems to operate in 
regimes where virial shock heating is expected to be unim- 
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portant, it seems to be a weak effect compared to the halo 
quenching which seems to have eventually overtaken M*- 
dependent quenching by z ~ 0. 



the weak correlation between M* and Mh at the high end 
for centrals and no correlation at all for satellites), which is 
the opposite of the behaviour of the quenched fraction (see 
Fig. Illl and the discussion in £|5.3|) . 



7 SUMMARY AND CONCLUSIONS 

Using the large sample size of the SDSS at z ~ and the 
multi-wavelength nature of the AEGIS dataset at z ~ 1 
we have examined the relations between stellar mass, halo 
mass, environment density and distance of satellites from 
the group centre over the redshift range < z < 1. We 
showed that there are two keys to understanding these rela- 
tions. First, the density estimator to the JVth-nearest neigh- 
bour behaves in two modes, one in which the number of 
observed group members is less than N so that the den- 
sity measures the distance to the nearest haloes, and the 
other in which the number is greater than N so that the 
density measures distances within a given halo. These are 
the cross-halo and single-halo modes respectively, and they 
are typically <5jv <C 1 and Sn 3> 1 for a given M*. Second, 
the relations between M» , Mh and Sn behave differently for 
satellites than for central galaxies. We argued that for satel- 
lites, the projected distance from the halo centre D, which 
anti-correlates with Sn in the single-halo mode, is a better 
measure of "environment" than 5n because the former does 
not behave in two modes that depend on both A and the 
limiting magnitude of the survey. D measures density, but 
also time spent in the cluster. 

Armed with the understanding of these relations, we 
explored quenching as a function of the quantities Mh, M„, 
Sn, and D. We used the wealth of SFR data from the SDSS 
and the AEGIS surveys to measure quenching as the fraction 
of low-SFR galaxies, rather than the fraction of red galaxies 
which can be severely affected by dust. We separated the 
analysis into samples of central and satellite galaxies since 
the predictions for the behaviour of the quenched fraction 
differ for these populations. 

We also discussed the effects of the halo mass function 
at z ~ 1 having fewer massive haloes compared to z ~ on 
the relations between M*, Mh and Sn, and on the behaviour 
of the quenched fraction as a function of these quantities. 

Our main results are listed in detail below: 

(i) Colour: an imperfect surrogate for quenching 

About one third (30% - volume weighted) of the red se- 
quence of the colour-magnitude diagram for the SDSS lies 
on the SF sequence. These dusty red star forming galaxies 
tend to be the most massive star forming galaxies and are 
on average 0.4 dex more massive than the whole population 
of SF sequence galaxies (see Fig. [3] and the discussion in 
jpS. T hese results are consistent with those of iMaller et al.l 
(|2009r ) who find that a third of red sequence galaxies move 
to the blue cloud after an inclination-dependent dust correc- 
tion, and that this dust correction is larger for galaxies with 
brighter Mk- The presence of dust-reddened galaxies on the 
red sequence means that colour does not correlate perfectly 
with star-formation quenching, which we showed has impor- 
tant implications for quenching correlations vs. stellar and 
halo mass. For both centrals and satellites, the red fraction 
will appear to correlate more with M* (due to the massive- 
ness of dusty star formers) and not much with Mh (due to 



(ii) Relations between Mh, M,, Sn and D 

(a) Environment density as measured by Sn, i.e., the 
density based on the distance to the Ath nearest neigh- 
bour, is characterised by two modes of behaviour, a cross- 
halo mode in which the number of observed group mem- 
bers is fewer than N and a single-halo mode where the 
membership is greater than A. 

(b) For central galaxies at z ~ 0, M* correlates with 
Mh below about M h ~ few x 10 12 M Q above which M* 
o nly weakly depend s on Mh, consistent with the findings 
of lYanget all J2009I ) (Fig.|IJi). At z ~ 1, the relation be- 
tween M, and Mh coincides with the relation at z ~ 
at Mh ~ 1O 12 M0, but rises above the z ~ relation for 
haloes more massive than that (Fig. [12b). This is con- 
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sistent with findings of lLeauthaud et al.l 
halo growth around central galaxies which grow very little 
in A/*. Most (85%) of all galaxies at z ~ 1 reside in haloes 
smaller than M h = 10 12 ' 5 M Q . 

(c) Central galaxies in the single-halo mode at z ~ 
are characterised by high <5g and high Mh, while those in 
the cross-halo mode exhibit little relation between £5 and 
Mh (Fig. \4%>). Central galaxies at z ~ 1 behave similarly, 
but the number of observed group members rarely exceeds 
A = 3 (Fig. Hit). 

(d) For satellite galaxies at 2 ~ in the single-halo 
mode (which is 58% of these satellites), M* does not 
strongly correlate with Mh- In the cross-halo mode for 
satellites, M* correlates with Mh as expected from statis- 
tical and dynamical friction constraints (Fig. 0>,c). 

(e) For satellite galaxies in the single-halo mode at 
z ~ 0, 5$ correlates with Mh (Fig. \5i>), as expected from 
considerations of the universal halo profile and from the 
predictions of the spherical collapse model of a virialised 
halo combined with the known halo occupation function. 

(f) For satellites at z ~ 0, 85 anti-correlates with 
D = d pro j/_Rvi r in the single-halo mode but shows no de- 
pendence in the cross-halo mode (Fig. [6|. This is consis- 
tent with the expected statistical behaviour of 85. 

(iii) Quenching, Mass, Density and Distance 

(a) For central galaxies at z ~ 0, the quenched fraction 
appears more strongly correlated with Mh (at fixed M*) 
than with M* (at fixed Mh) for centrals more massive 
than M, ~ 1O 1O ' 4 M (Fig. lib). 

(b) At z ~ 1, where most of the galaxies reside in 
haloes less massive than the halo-quenching scale, quench- 
ing correlates with M* rather than Mh for haloes below 
~ 1O 12 ' 5 M (which host 85% of these centrals), but the 
quenched fraction remains low ( < 0.2) (Fig. 113b ). This 
behaviour may signal a M«-dependent quenching mode at 
high redshift that is overtaken by Mh-dependent quench- 
ing as the haloes grow by 2 ~ 0. 

(c) For central galaxies at 2 ~ 0, in the <5s-M* plane, 
the quenched fraction increases with both M* and 85 (es- 
pecially at high M„) in such a way that is consistent with 
the increase in mean Mh with both of these quantities 
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(compare Fig. [8ji with Fig. [4}:). In other words, the in- 
crease of quenching with both M* and 85 for centrals 
largely reflects the increase of quenching with Mh alone. 
At z ~ 1 the central galaxies behave in a similar way 
(compare Fig. [T3b with Fig. [I2b 1 . 

(d) For satellite galaxies at z ~ 0, the quenched frac- 
tion increases with Mh and decreases with the distance 
to the group centre D (Fig. llOb ). The combination of Mh 
and D is an excellent predictor of quenching for satellites, 
producing quenched fractions as high as 95%. Mh and D 
are preferred over M, and <5g as predictors of quenching 
because 1) 85 behaves differently in two modes and suf- 
fers from redshift effects making quenching trends with 
85 difficult to interpret, and 2) quenching follows average 
Mh better than average M„ (compare Fig. [SJ; with Fig. [5£ 
and Fig. [9j; with Fig. \5fc). M* -dependent quenching seems 
to matter more than Mh only at large distances from the 
group centre (where satellites are likely to have just re- 
cently fallen in) and may reflect a quenching dependence 
on the mass of the sub-halo (compare Fig. 1 10b with b). But 
nearer to the group centre than ~ 0.2_R v i r , the quenching 
depends strongly on host Mh and almost not at all on 
satellite M». 

Our main conclusion is that the quenched fraction of 
galaxies at z = — 1 seems to be primarily driven by 
the dark-matter halo mass. This is consistent with the 
theoretical prediction for an accretion shutdown due to 
virial shock heating once the halo grows above a thresh- 
old mass, which possibly triggers other quenching mecha- 
nism s , such as stellar and AGN feedback (jBirnboim fe Dekell 
120031 : iDekel fe Birnboim] bOOrJ ). For satellite galaxies, the 
quenching also depends on the proximity to the halo cen- 
ter, which reflects the local density and the history of the 
satellite within the halo. Stellar mass is only a secondary fac- 
tor, and Sn is a poor measure of environment. As discussed 
in £15.31 o ur results are consisten t with those of previous 
work (eg.. rWeinmann et al.ll20 06a; Bla nton fc Roweisll2007l ; 
lHansen et aLll2009l ; iKimm et aLll2009l ; IWetzel et al.ll2012ah . 
Our re sults can be brought into agreement with the conclu- 
sion of IPeng et al.l (|2010l . 120111 ) once stellar mass is substi- 
tuted by halo mass. The difference is driven by our proper 
use of SFR rather than color to define quenching, and it has 
the advantage of unifying all phenomena under a common 
halo quenching model. However, marginal evidence is also 
found for an M* -related quenching mode at earlier epochs 
when Mh is low. It will be valuable to study the quenching 
mechanism with comparisons to galaxy properties beyond 
M», such as bulge-to-disk ratio and central density (Cheung 
et al., Fang et al. and Barro et al., in preparation). Combined 
with the accretion shutdown associated with halo mass, in- 
ternal mechanisms related to the bulge and the central black 
hole ma y serve as the direct cause for shutdo wn in star for- 
mation l|Dekel et al.ll2009l ; iMartig et alj|2009h . 
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